Interaction effects in 2D electron gas in a random magnetic field: Implications for 
composite fermions and quantum critical point 
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We consider a clean two-dimensional interacting electron gas subject to a random perpendicular 
magnetic field, h(r). The field is nonquantizing, in the sense, that A/^-a typical flux into the area A F 
in the units of the flux quantum (Ap is the de Broglie wavelength) is small, Nh -C 1. If the spacial 
scale, £, of change of h(r) is much larger than A F , the electrons move along semiclassical trajectories. 
We demonstrate that a weak field-induced curving of the trajectories affects the interaction-induced 
electron lifetime in a singular fashion: it gives rise to the correction to the lifetime with a very 
sharp energy dependence. The correction persists within the interval cu ~ ujq — E^Nl' 3 much 
smaller than the Fermi energy, E F . It emerges in the third order in the interaction strength; the 
underlying physics is that a small phase volume ~ (uj/Ep) 1 ^ 2 for scattering processes, involving 
two electron-hole pairs, is suppressed by curving. Even more surprising effect that we find is that 
disorder-averaged interaction correction to the density of states, 5v(lo), exhibits oscillatory behavior, 

periodic in (u) /uo^ 3 ^ 2 . In our calculations of interaction corrections random field is incorporated 
via the phases of the Green functions in the coordinate space. We discuss the relevance of the new 
low-energy scale for realizations of a smooth random field in composite fermions and in disordered 
phase of spin-fermion model of ferromagnetic quantum criticality. 

PACS numbers: 71.10.Pm, 71. 10. Ay, 71.70.Di, 73.40.Gk, 73.43.Nq 



I. INTRODUCTION 

Electron-electron interactions are strongly modified 
when electrons move diffusively^. Resulting enhance- 
ment of the interactions leads, in two dimensions, to a 
divergent correction to the density of states^, Sv(ui). 
When electrons move ballistically and are scattered by 
point impurities, the anomaly persists, although it has a 
different underlying scenario^. 

Within this scenario, individual impurities (unlike the 
diffusive casei^) are responsible for the ballistic zero- 
bias anomaly by virtue of the following process. Static 
screening of each impurity by the Fermi sea creates a 
Friedel oscillation of the electron density with a period, 
A F /2, where A F is the de Broglie wavelength. Then 
the amplitude of combined scattering from the impu- 
rity and the Friedel oscillation, which it created, exhibits 
anomalous behavior— when the scattering angle is cither 
or 7T. Energy, ui, of the scattered electron, measured 
from the Fermi level, E F , defines the angular interval, 

1/2 

~ (ui/E F ) , within which the scattering is enhanced. 
This enhancement translates into 5v{oj) oc lnw correc- 
tion to the density of states. 

In a diagrammatic language, creation of the Friedel 
oscillation is described by a static polarization bubble. 
We note in passing, that the same polarization bubble 
at finite frequency, w, is responsible for the lifetime of 
electron of energy ~ u> with respect to creation of an 
electron- hole pair. 

It is knowns that, in perfectly clean electron gas, finite- 
range interactions do not cause any anomaly in 8v{uj). 
Then a natural question to ask is whether or not the 
anomalous behavior of 5v(u>) holds when a weak disor- 
der is not point-like, as in Ref. H but, instead, smooth. 



Finding an answer to this question is the main objec- 
tive of the present paper. For concreteness we choose a 
particular case of 2D electron gas in a smooth random 
magnetic field, although our main results apply to the 
arbitrary smooth disorder. 

Historically, the interest to the problem of 2D elec- 
tron motion in a random static magnetic field first 
emerged in connection with a gauge field description 
of the correlated spin systems^i Later this interest 
was stimulated by the notion that electron density vari- 
ations near the half-filling of the lowest Landau level 
reduces to random magnetic field acting on compos- 
ite fermions^. Another motivation was the possibility 
to realize an inhomogeneous magnetic field, acting on 
2D electrons, artificiallyi^^H^ds^OLIs. For non- 
interacting electrons, this motion has been studied the- 
oretically in Refs.tTliiilllii^^ 

|32| . In the present paper we trace how the perturbation of 
electron motion by a smooth random field affects the in- 
teraction corrections to the single-particle characteristics 
of the electron gas. 

In Rcfs. 5.6.7 the averaging over static random field 
was carried out with the help of the path integral ap- 
proach originally employed for diffusively moving elec- 
trons in a noisy environment (see also Refs. HHH!). A 
crucial fact that ensures the effectiveness of this approach 
is that the field is assumed to be 5-correlated. In fact, the 
correlation radius must be even smaller than A F . How- 
ever, in r ealizations^2J£JJJ5ii2Jiii5iiSiiLi£ mentioned 
above, the spatial scale of change of the random field 
in much bigger than A F . This leads to a completely 
different, semiclassical, picture of the electron motion, 
when only the paths close to the classical trajectories 
are relevant. In the present paper we consider only this 
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limit. Semiclassical character of motion suggests the way 
in which to perform the averaging over disorder realiza- 
tions. Namely, the equation of motion can be first solved 
for a given realization, while averaging over realizations 
is carried out at the last step. This order is opposite to 
Refs. where averaging was carried out in the gen- 

eral expression for the Green function after it was cast in 
the form of a path integral. 

It might seem counterintuitive that any smooth dis- 
order could generate a low-frequency scale for the inter- 
action effects. Indeed, smooth random field (including 
magnetic) does not produce Friedel oscillations, which 
are required for the anomaly 3 - to develop. In a formal 
language, there are no static bubbles in the diagrams for 
the interaction correction to the self-energy. More pre- 
cisely, in the smooth random field, they are exponentially 
suppressed. We will, however, demonstrate that the low- 
frequency scale emerges from dynamic bubbles after they 
are modified by a smooth disorder. 

The new low-to scale shows up in the virtual processes 
involving more than one electron- hole pair, i.e., two or 
more bubbles. This is because the momenta of states, 
involved in these processes, are strongly correlated, as 
was first pointed out in Ref. l3rl Namely, these momenta 
are either almost parallel or almost antiparallel to each 
other. It is this correlation in momenta directions that 
is affected by the smooth random magnetic field. By 
suppressing the correlation, random field gives rise to the 
low-w feature in 5v(uj). Clearly, both the height and the 
width of the feature, depend on the magnitude of the 
random field. The above argument makes it clear why 
the low-w scale does not emerge on the level of a single 
bubble, modified by the random field. The reason is that 
the single bubble describes excitation of a single pair; 
there is no strong restriction on the momenta directions 
in this process. 

Once the mechanism of nontrivial interplay of smooth 
disorder and interactions is identified, the following ques- 
tions arise: what is the shape of the anomaly in 5i>(lu), 
and how it depends on the strength and the correlation 
radius of the random field? To address these questions 
we develop a systematic approach to the calculation of 
interaction corrections in a smooth random field. The 
key element of our approach is incorporating the ac- 
tion along the curved semiclassical trajectories into the 
phases of the Green functions. Our calculation reveal 
a surprising fact, which could not be expected on the 
basis of the above qualitative consideration. It turns 
out that disorder-averaged correction, (<5^(tt>)), exhibits 
an oscillatory behavior. Oscillations emerge when two 
pairs, participating in one of the possible processes giv- 
ing rise to <5z/, are strongly correlated with each other. 
As an example consider the process, involving creation of 
the electron-hole pair, rescattering within the pair, and 
its subsequent annihilation. In this process, oscillations 
come from electron-electron scattering events that hap- 
pen at the points, located on a straight line and at equal 
distances. To the best of our knowledge, this is the first 



example when disorder does not suppress, but on the 
contrary, brings about the oscillations. 

Therefore, as we demonstrate in the present paper, 
anomaly in the density of states is created by smooth 
spatial variation of the magnetic field, even though this 
variation does not produce Friedel oscillations. Although 
modification of the Friedel oscillations from a point-like 
impurity by a smooth random field is not directly related 
to our situation with no impurities, this problem is still 
useful for gaining a qualitative understanding. Indeed, 
the relevant random-field-induced length scales, in our 
clean case, emerge in this problem as well. For this rea- 
son we start with the study of suppression of the Friedel 
oscillations by the random field, prior to the analysis of 
the interaction corrections in the random field. 

We are not aware of literature on disorder-induced 
smearing of the Friedel oscillations 3 ^. However, a closely 
related issue of smearing of Rudcrman-Kittel-Kasuya- 
Yosida (RKKY) interaction between the localized spins 
by the disorder, has a long history^ 3 ^^! 4 ^ 3 .^. It ^ 
easy to see2& that a short-range disorder suppresses expo- 
nentially the average RKKY interaction. Howeve r 39 i 40 ' 41 , 
the average interaction does not represent the actual 
value of exchange in a given realization. This is due to the 
fast oscillations of the exchange with distance. The typ- 
ical magnitude of the exchange can be inferred from the 
averaging of the square of the RKKY interactio n 39 ' 40 ' 41 ; 
this average is suppressed by the disorder only as a power 
law. 

In this paper we demonstrate that the decay of the 
averaged Friedel oscillations in the presence of a smooth 
disorder is quite nontrivial. In particular, when the field 
is strong enough, the average, in contrast to Ref. |3g, falls 
off with distance as a power law. We would like to note 
that recently the notion of averaged Friedel oscillations 
became meaningful. This is because the possibility of 
visualization of a single-impurity-induced oscillation had 
been demonstrated experimentallyi 5 .^^^^ 5 ^^^. 
The role of averaging can be then played by slow tem- 
poral fluctuations of the environment. Since experimen- 
tal advancesi^^^^iSi^iS 3 . were reported for cor- 
related systems, recent theoretical studie s 54 ' 55 ' 56 ' 57 ad- 
dressed the Friedel oscillations created by a single impu- 
rity in such systems. 

The paper is organized as follows. In Section [II] pos- 
sible regimes of electron motion in a random magnetic 
field are identified. In Section IlIII we summarize our re- 
sults on Friedel oscillations and interaction correction to 
the density of states for weak random field, i.e., for the 
field, in which the straight-line electron trajectories are 
weakly perturbed by the field. Subsequent Sections IIVI 
HXI are devoted to the derivation of the results, outlined in 
Section HTH Finally, In Section [X] we translate our results 
into predictions for experimentally observable quantities 
in two prominent situations: composite fermions in half- 
field Landau level and electrons interacting with critical 
magnetic fluctuations near quantum critical point. De- 
tails of some of the calculations are presented in Appen- 



I 



II 



e 




III 



© 









r 





V : 




u=k F R L 



FIG. 2: (Color online) Parametric regions for the regimes I, 



II, and III. Dashed line, 



-1/3 



, separates slow and fast 



power-law decays of the averaged Friedel oscillations within 
the regime I: the oscillations fall off as 1/r 2 to the left from 
the dashed line, and as 1/r 7 / 2 to the right from the dashed 
line. 



FIG. 1: (Color online) Types of semiclassical trajectories be- 
tween two points separated by a distance, r, in a random 
magnetic field: in the regime I the trajectories are of "arc"- 
type; in the regime II the trajectories are of "snake"-type; 
regime III corresponds to a drifting Larmour circle. 



dices A-F. 



II. REGIMES OF ELECTRON MOTION IN A 
RANDOM MAGNETIC FIELD 

Let r = (a;, y) be the coordinates of the 2D electron. 
Random magnetic field along z-direction is characterized 
by the correlator 

{h{v)h{v')) = hlK{\v-v'\/i), K(0) = 1, (2.1) 

where ho is the r.m.s magnetic field and £ is the correla- 
tion radius. Throughout the paper we will assume that 
the random field is slow- fluctuating, in the sense, that £ 
is much bigger than the de Broglie wavelength A F , the 
case opposite to the limit £ — > considered in Refs.^HQ- 
In terms of semiclassical description, different regimes of 
motion are classified according to the classical electron 
trajectory, which begins at the origin and ends at point 
r. One should distinguish three different regimes, as il- 
lustrated in Fig. 1. 

(z) short-distance regime (regime I in Fig. 1). The tra- 
jectory is of the arc-type. For this regime to realize, two 
conditions must be met. Firstly, the change of magnetic 
field over the distance, r, should be negligible, i.e., r -C £. 
Secondly, the curving of electron trajectory in the locally 
constant magnetic field must be relative small. The mea- 
sure of this curving is r/i? L , where i? L = hck F /eho is 
the Larmour radius in the field, ho, and k F = 2tt/X f is 
the Fermi momentum. Thus the short-distance regime 
corresponds to r <C £, i? L - 



(ii) "weak-field" long-distance regime (regime II in 
Fig. 1). The trajectory is of the snafce-type. One con- 
dition for this regime is that magnetic field changes sign 
many times within the distance, r, i.e., r ^> £. The other 
is that within each interval of length ~ £ the curving of 
the trajectory is weak, i.e., £ <C i? L - 
(ra) " strong-field" long-distance regime (regime III in 
Fig. 1). Electron executes many full Larmour circles 
before arriving to the point r. The conditions for this 
regime are fl L <r and ii L <(. 

Note, that the last two regimes correspond to the 
"semiclassical" and "strong" random magnetic field 
regimes in the language of Ref. |26|. In order to accom- 
modate all three regimes within a single diagram, it is 
convenient to introduce the dimcnsionlcss parameters 



u •■ 
v ■ 



k F R L 



(chk F /eho) 
k F rJ\f h , 



(2.2) 



where Mh < 1 is the flux of the field ho into the area 
A F (in the unites of the flux quantum) . Then the regime 
I is defined by the lines u = fc F £ and v — k F £/u, see 
Fig. 2. The regime III is separated from the regime I 
by the line v = 1, and from the regime II by the line 
u = fc F £. Finally, the dashed region u < 1 in Fig. 2 
corresponds to quantizing magnetic field. The diagram 
Fig. 2 is compiled for fc F £ 3> 1, so it does not reflect 
white-noise regime, fc F £ <§; 1, of Refs. [^laB 



III. MAIN RESULTS 

A. Friedel oscillations 

The simplest manifestation of the interplay of external 
field and electron-electron interactions shows up in spa- 
tial response of the electron gas to a point-like impurity, 
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or, in other words, in Friedel oscillations. Denote with 
Ui mp (r) the short-range potential of the impurity. In the 
presence of interaction, V(r — ri), the effective electro- 
static potential in a clean electron gas falls off with r as 
V H (r) oc sin(2fc F r)/r 2 in a zero field. In Rcfj££ we had 
demonstrated that in a constant magnetic field, h = ho, 
this behavior modifies to 



V B (r) = - 



iy gV(2k F ) 
2nr 2 



sin 



2n/p 



(Por) 3 



12 



(3.1) 



where the characteristic momentum, po, is defined as 

/ \ 2//3 

k F I ho } 

Po 



(fc F i? L ) 2 / 3 



kl'^o, 



(3.2) 



IS 



where $o is the flux quantum. In Eq. (|3.ip vo = m/nh 
the free electron density of states, V(2k F ) is the Fourier 
component of V(r), and the parameter g is defined as 
g = J Ui mp (r) dr. Eq. (|3.5p is valid within the domain 
fcjJT 1 < r < i? L , so that (por) 3 /12 in the argument of sine 
does not exceed the the main term, 2k F r. As follows from 
Eq. (|3.2| . the characteristic length scale, 



/3 / $o 



2/3 



(3.3) 



ri = — = />v I , 

Po \ ho . 

defied by po, is intermediate between i? L and A F , so that 

i? L >r I >l/fc F . (3.4) 

We see from Eq. (13. ip that only the phase of the Friedel 
oscillations is affected by the constant field, while the 
magnitude still falls off as 1/r 2 . The randomness of 
h(x, y) results in randomness of the field-induced phase 
of the oscillations. This, in turn, translates into a faster 
decay of disorder-averaged oscillations. To quantify the 
behavior of the average (V s (r)), we rewrite it the form 



(V H (r)) = - 



v gV{2k F ) 
2irr 2 



F(r) sin 



2k F 



4>{r) 



(3.5) 



so that F(r) describes the decay of the magnitude of the 
disorder-averaged oscillations. For a given distance, r, 
the character of the phase randomization is different in 
the regimes I and II. In regime I, we have £ S> r, and 
thus the relevant scale for the decay of (V B (r)) is r : . In 
Section [V] we find that in this regime the magnitude, 
F 1 , and the phase, <j> u are the following functions of the 
dimensionless ratio x — r/fi 



1 



(1 + x 6 ) 174 ' 
4>i(x) = — arctan 



(3.6) 



In regime II, with snake-like trajectories, Fig. [TJ the 
sign of random field changes many, ~ r/£ 3> 1, times 
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FIG. 3: (Color online) Friedel Oscillations of the potential 
created by an impurity located at the origin. Thick line: aver- 
aged Friedel oscillations in regime I is plotted from Eqs. (|3.5|l , 
(|3.6p . Thin line: oscillations in the absence of the random field 
[Eq. (33J with p = 0]. 



within the distance, r. As demonstrated in Section |V1 in 
this regime we have 



F u (x) 



V2a 



1/2 



(jhi{x) 



1 + |x 4 v cosh 2 x — cos 2 x 



arctan 
arctan 



(3.7) 



1 — cot x tanh x 



(3.8) 



where x — r/r u , with r„ defined as 



r u = V 



, v 1/2 

U V 



(3.9) 



In Eq. (|3.9|) the numerical factor, rj, depends on the func- 
tional form of the correlator Eq. (|2.ip and will be defined 
in Section [V] 

In conclusion of this subsection we point out that the 
actual character of the decay of Friedel oscillations with 
distance is governed by the following dimensionless com- 
bination of parameters, ho, and, £, in the correlator 
Eq. (|2.ip of the random field 



2e3 



<S> 2 k F 



(3.10) 



For e 3> 1, i.e., for strong random field, the averaged 
oscillations decay with r according to Eq. (|3.6p in the 
regime I. This is because for £> 1 we have pot; 3> 1. In 
the opposite limit of a weak random field, e <ti 1, we have 
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Pot, <C 1, so that the scale Pq 1 is irrelevant, and also no 
dephasing takes place within the distance, £. Thus, the 
characteristic decay length, r„ ~ ^/e 1 / 2 , is much larger 
than £. This automatically guarantees that k F r u S> 1. 



B. Tunnel density of states 

Two spatial scales, r : and r„, define two energy scales, 

2/3 



Wl = 







' /10 \ 






Til 


-£p 


(e /2 h Q 



E F Af, 



2/3 



(3.11) 



E F {k^) 1/2 N h 



As shown below, these scales manifest themselves in 
the anomalous behavior of the density of states in the 
third order in the electron-electron interaction parame- 
ter, vqV . More specifically, in the regime I, the bare 
density of states, vq, acquires a correction 8pi(u>) ~ 
^(^o^) 3 ( w o/ E F ) 3 ^ 2 X (uj /ujq) , where E F is the Fermi en- 
ergy. In the regime II the correction has a similar form 
5u u (u}) ~ M"oV) 3 {lui/E f ) 3/2 J {lu/lu!). Both func- 
tions, T{z) and J{z), have characteristic magnitude and 
scale ~ 1. Moreover, they exhibit quite a "lively" be- 
havior. In particular, a zero-bias anomaly, 8v 1 {uj), falls 
off at w > wq with aperiodic oscillations, i.e., T{z) has 
a contribution oc sin (2 8 / 3 V3z)z~ 3 / 4 exp {-2 8 / 3 z} for 
2 > 1, The origin of the oscillations is the power-law 
decay of F 1 (x), given by Eq. (|3.6p . and the brunch-point, 
x = e™/ 6 . 

The contribution, 8v n (uS), also has a non-monotonic 
behavior, despite the fact that F u (x) falls off exponen- 
tially, as exp(— r/m) [see Eq. (I3.7|) ]. 

It is instructive to trace the evolution of the zero-bias 
anomaly upon increasing the magnitude of the random 
field, Hq. This evolution is governed by parameter, e, 
Eq. (|3.10p . While e remains smaller than one, where 
the regime II applies, the anomaly is described by the 
function J(oj/ujx) and broadens with ho as v F /r u (ho) oc 
ho. Upon further increasing ho, when e exceeds one, the 
crossover to the regime I takes place. Zero-bias anomaly 
is then described by T{u>/u>o)', it broadens with ho as 

Vp/r^ho) oc h^ 3 , and develops oscillations. The fact that 
oscillations in 8v(u>) emerge upon strengthening disorder 
might seem counterintuitive. This issue will be discussed 
in details in Section IVTil 

In a zero magnetic field, an intimate relation be- 
tween impurity-induced Friedel oscillations and the zero- 
bias anomaly was first established in Ref. H|. Namely, 
it was demonstrated that for short-range interaction 
5v(lo)/vo ~ (voV/E f t) lnw, where 1/r = fo7r<? fMmp is 
the electron scattering rate by the impurities, and n- lmp 
is the impurity concentration. This anomaly is of the 
first order in V. A non-trivial question is whether or not 
the modification, Eq. (|3.1|) . in a constant magnetic field 
results in field dependence of the density of states in this 
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FIG. 4: (Color online) Magnetic-field-induced contribution 
Eq. (|7.25[) to the ballistic zero-bias anomaly Eq. (|3.12[) . Field- 
dependent correction, [Sv(oj, h) — 5v(ui, 0)] /vq, in the units 
(foF) 2 (wo/S F ) 1/2 (£ , F r)- 1 is plotted versus dimensionless en- 



ergy 2 2/3 (w/wo), where w = (2E F ) 1/3 cu 2 c /3 
the cyclotron frequency. 



» uj c , and uj c 



order. In other words, whether or not a weak magnetic 
field introduces a cutoff of In a; at small u>. The answer 
to this question is negative. In Ref. 58 it was demon- 
strated that sensitivity of 8v{u>) to a weak magnetic field 
indeed emerges, but in the second order in vqV (how- 
ever, still in the first order in 1/r). The field-dependent 
correction, [8v(u>, h) — 8v(lo, 0)] jvo, has a characteristic 
frequency scale, to = loq. It is interesting to note that, at 
lj ^> loo, this impurity-induced correction has an oscillat- 
ing character 



8v(cj, h) - 8v{uj, 0) (voV) 2 ( tu \ 1/2 



vo 



E f t V E F J \uo 



(3.12) 



The dimensionless function, P, has the following large- a; 
asymptote 



P(x) oc 



1 



,3/4 



r 3/2 



3V3 



(3.13) 



In Fig. 2] we show the oscillating correction to the density 
of states; the form of the function P (ui/ujq) is addressed 
in Section lVIII Technically, the derivation of Eqs. (|3.12l) . 
(13.13)) is quite analogous to the derivation of the oscilla- 
tory 8v in the random field in the regime I. For this reason 
we will outline this derivation in Section [VIII 



IV. POLARIZATION OPERATOR IN A 
RANDOM MAGNETIC FIELD 

Friedel oscillations, V H (r), created by a point-like im- 
purity, and the ballistic zero-bias anomaly originating 
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from these oscillations^ are intimately related to the 
Kohn anomaly in the polarization operator, n(q), of a 
clean electron gas near q = 2k F . In two dimensions, this 
anomaly behaves as~L {q— 2fc F ) 1 / 2 , which translates into 
1/r 2 decay of the Friedel oscillations and a lnw correc- 
tion to the density of states. Suppression of the Friedel 
oscillations, V H (r), in a random field is a result of smear- 
ing of the Kohn anomaly in the momentum space. How- 
ever, since the momentum is not a good quantum number 
in the presence of the random field, it is much more con- 
venient to study the field- induced suppression of V n (r) 
directly in the coordinate space. 



A. Evaluation in the coordinate space 

Polarization operator, nn(r,r'), is defined in a stan- 
dard way as 

H(r, r', n) = -i J ^G n > (r, r')Gn-n' (r', r). (4.1) 

Here Gn(r, r') denotes causal Green function, which coin- 
cides with the retarded, Gq(t,t'), or advanced G^(r,r') 
Green functions for £1 > and il < 0, respectively. At 
distances |r— r'| 3> fcpT 1 the polarization operator in coor- 
dinate space represents the sum ELj(r, w) and n 2 fc F (r, u>) 
of slow and rapidly oscillating parts 



n (^) = -^- M exp J- 



v Q h 3 . \ . f 2irrT\ 



2r 2 



xexpj^j. (4.3) 



Subindices and 2fc F emphasize that these parts come 
from small momenta and momenta close to 2fc F in TL(q), 
respectively. Eq. (|4.2[) emerges if one of the Green 
functions in Eq. (|4.1[) is retarded and the other is ad- 
vanced. Eq. (|4.3p corresponds to the case when the 
Green functions in Eq. (|4.ip are both advanced or both 
retarded^^. Derivation of Eqs. (|4.2|) . (|4.3|) is presented 
in Appendix A. In Eq. (|4. 3[) the function, 



= - 



sinhx 

in n2fc F describes the temperature damping. 



(4.4) 



Classical trajectory of an electron in a weak magnetic 
field is curved due to the Larmour motion even at the 
spatial scales much smaller than R L . As a result of this 
curving, the electron propagator, G(ri,r 2 ), between the 
points i"i and r 2 contains, in the semiclassical limit, a 
phase, k F C, where C is the length of the arc of a circle 
with the radius i? L , that connects the points ri and r 2 , 
see Fig. [5£i. Since the Friedel oscillations are related to 
the propagation from ri to r 2 and back, it is important 
that two arcs, corresponding to the opposite directions of 
propagation, define a finite area, A, so that the product 
G(ri, r 2 )G(r 2 , ri) should be multiplied by the Aharonov- 
Bohm phase factor, exp[ihoA/$o]. Then the phase, of 
this product is equal to 

2fc F r + 0(r) =2fc F £- h " A ^ r *\ (45) 
Simple geometrical relations, see Fig. [5^, yield 



r = \rx 
£ = RJ, 



-r 2 | 
A-- 



= 2i? L sin(<5/2), 
2i? 2 (<5-sin<5). 



Using this relation and assuming r « fi^ we find 



4>(r) 



12/c F $ 2 



(Po^) £ 
12 



(4.6) 



(4.7) 



At this point, we would like to note, that the conventional 

(4.2) 

way£2- of incorporating magnetic field into the semiclas- 
sical zero-field Green's function amounts to multiplying 
it by exp [(l/$o) / a • dl] , where the phase factor is the 
integral of the vector potential, a, along the straight line, 
connecting the points ri and r 2 . Such an incorpora- 
tion neglects the field-induced curvature of the electron 
trajectories, and thus does not capture the modification 
Eq. (|3.5[) of the Friedel oscillations in magnetic field. In- 
deed, the magnetic phase factors, introduced following 
Ref. [62j cancel out in the polarization operator. 

With phase, 4>(r), given by Eq. (|4.7[) . Friedel oscil- 
lations in a constant magnetic field acquire the form 5 ^ 
Eq. (|3.1[) . To sec this, we notice that, with accu- 
racy of a factor, g/2iT, the potential, V H (r) coincides 
with II 2fcF (r, 0). Then the additional phase Eq. fl~7|) 



transforms sin(2fc F r) into sin [2fc F r — (por) 3 /l2^ , as in 
Eq. (|3~TD . 

In Appendix B we present a rigorous derivation of 
Eq. (|3.1|) starting from exact electronic states in a con- 
stant magnetic field, as in Ref. [&3l 



B. Qualitative derivation for the constant field 

For a constant magnetic field, h(x,y) = /to, the phase, 
4>(r), in the argument of Eq. (|3 . 5[) can be inferred from 
the following simple qualitative consideration. 



C. Field-induced phase of the Green function: 
Analytical derivation in a spatially-inhomogeneous 
field 



Additional semiclassical phase, 5ipo^ r , of the Green 
function due to the random magnetic field, h(x,y), is 
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R i 



r 2 



R, 




a) 



b) 



FIG. 5: (a) Origin of the net "magnetic" phase Eq. (14. 7ft : two 
arcs, corresponding to the opposite directions of propagation, 
define a finite area, A. Aharonov-Bohm flux through this 
area makes the net phase negative; (b) Schematic illustration 
of the scattering processes giving rise to the additional phases 
Eqs. (EU), in the product Eq. (J775J). 



given by the following generalization of Eq. (|4.7 



<fy>0->r = 



k f r ( dy \ 1 f r 



h(x,0), 



where the first term comes from the elongation of the 
trajectory in magnetic field. The second term describes 
the Aharonov-Bohm flux into the area restricted by the 
curve y(x) and the z-axis. In Eq. (|4.8p we assumed that 
the field does not change along the y-axis. This is the 
case when the maximal y is smaller than the correlation 
radius, £, of the random field. The condition y < £ is 
met in the regime of the "arcs" and the regime of the 
"snakes", see Figs.[Q[l 

In Eq. (|4.8|) we have also assumed that the magnitude 
of the de Broglie wavelength of the electron does not 
change along the trajectory. This can be justified from 
the equations of motion 



d2 y e hi n 



m—7r = h(x,0) — . 

dt 2 c y " 



dx 
It 
dy 
dt ' 



(4.9) 



It follows from Eq. (14. 9|) that the energy of electron 
Y [(dx/dt) 2 + (dy/dt) 2 ] is conserved even if magnetic field 
changes with coordinates. 

The most important step that allows to find dipo— >r 
analytically, is that in the regimes I and II in Fig. 1 we 
can replace dx/dt by v F and set t = x/v F in the rhs of 
Eq. P~!?)) . This allows to replace d 2 y/dt 2 by vjd 2 y/dx 2 . 
Then the first of the equations yields 



m v 



2 d y 

dx 2 



ev T . 



h(x,0). 



(4.10) 



Integrating this equation, we obtain 



dy _ 

dx mcv 



dx'h{x\ 0) + C. 



(4.11) 



F JO 



The constant, C, should be found from the conditions: 
y(0) = 0, and y(r) — 0, leading to 



C = 



dx' / dx"h(x",0) 
mcv F r J J 



dx'A(x'), (4.12) 



mcv F r jo 

where we have introduced an auxiliary function 



A(ar) = / dx'h(x',0). 
Jo 



(4.13) 



The meaning of A(x) is the y-projection of the vector 
potential. Substituting Eq. (|4.12|) back into Eq. (14. lip , 
we find 



dy _ e 
dx mcv F 



dx'h(x',0) 



dx' / dx"h(x",0) 
ii Jo 



mcv F 



1 



A(x) / dxA(x) 



(4.14) 



With the help of Eq. (|4.14[) one can express the first term 
in additional phase Eq. (|4.8[) in terms of A(a;). It turns 
out that the second term in Eq. (14. 8|) exceeds twice the 
first term. To see this, one should multiply the first of 
equations Eq. (|4.9p by y{x) and integrate over x 



dx h{x,0) y{x). (4.15) 



dx yi^x^j — — 

(i dx 2 $ fc F Jn 



The rhs of Eq. (|4.15p is the second term in Eq. 
The lhs of Eq. (|4.15p can be related to the first term in 
Eq. (|4.8p upon integration by parts 



dx y{x) 



dx 2 



dx 



Finally, we get 



K fdx( d ^ 
2 Jo \dx 



(4.16) 



(4.17) 



-1-1 / dxA 2 (x)-- 

*0 fc F \ Jo r 



1 



dx A(x) 



It is convenient to rewrite the final result Eq. (I4.17P di- 
rectly in terms of the random field, h(r). Substituting 
Eq. (|4TT5|1 into Eq. (|4"T7|) . we obtain 
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6<Po-+, = J dr x J dr^h^nir^r^jhir^), (4.18) 

where the dimensionless kernel 7?.(r 1 ,r 2 ) is defined as 

^(ri,r 2 )=£<%i)% 2 ) (4.19) 

- — - - x 2 0(x 2 - x\) - x\ 6(x\ - x 2 ) 
r 

Note, that for the constant field h(x, y) = h , evaluation 
of Eq. (|4.18p using the kernel Eq. (|4.19[) reproduces the 
result Eq. (|4.7[) . as expected. 



V. DISORDER-SMEARED FRIEDEL 
OSCILLATIONS IN DIFFERENT REGIMES 

Smearing of the Friedel oscillations in the random field, 
h(r), originates from the randomness of the phase, ifo->r, 
which is related to h(v) via Eqs. (|4.18[) . (|4.19[) . Quan- 
titatively, the magnitude, F(r), and the phase, (f>(r), of 
smeared Friedel oscillations Eq. (|3.5p are determined by 
the following averages 



T x (r) = Im(e 2i5 ^> T 2 (r) = Re( e 2i ^— 



(5.1) 

Then F(r) and cf>(r) are related to the functions Ti(r) 
and T 2 (r) as 



F(r) = 



0(r) = arctan 



T a (r) 



(5.2) 



In this Section the averages Eq. (|5.ip will be calculated 
separately for the regime of "arcs" and the regime of 
" snakes" . 



A. Regime I 

In the regime of "arcs" we have r <C £, so that the field 
is almost constant within the interval (0, r) and is equal 
to its "local" value. For this reason, we can perform the 
averaging of exp [2i</3o— >r] over realizations of the random 
field, h(x,y), explicitly, without specifying the form of 
the correlator, K(r/£). This is because we can first set 
h(x,y) = const in exp [2iipo^ r ] , and then make use of 
the fact that the distribution function of the local field is 
Gaussian^. Characteristic spatial scale, r x , for F(r) and 
4>{r) immediately follows from Eq. (|4.18[) upon setting 
/i(0, r) = ho, and requiring 26(po^ r = 1. This yields 
r r = 2 2 / 3 3 1/ ' 3 /po, where po is given by Eq. 



1. Random magnetic field 

As discussed above, we start with Friedel oscillations 
in a constant local magnetic field, h, for which we know 
that 



F I (r,h) = l, 4>i(r,h) 



(5.3) 



where po = k F (ui c / E F ) 2 ' 3 , and uj c — eho/mc is the cy- 
clotron frequency in the field, ho. In Eq. (|5.3[) the pa- 
rameter, e r , is defined as 



h 2 r 3 
12$2fc p 



{Par) 3 
12 



(5.4) 



To find the form of the averaged Friedel oscillation in the 
regime I, in which po£ <1, we have to simply substitute 
the "local" value, h, of magnetic field into Eq. (|3.1[) . i.e., 
replace Pq by Pgh 2 jh\, and perform the gaussian averag- 
ing over the distribution of the local field. This averaging 
can be carried out analytically with the use of identity 



oo 



* dx 



• 2 , 

cos(e r 



! + /3) = Ti(e r ) cos (3 - T 2 (e r ) sin/3, 



(5.5) 



where the functions Ti and T 2 for this case assume the 
following forms 



(5.6) 





(5.7) 



Using Eq. (|5.2p . we recover from Eqs. (|5.6[) . (|5.7[) the final 
result Eq. f|3. 6|) for the magnitude F^r/r^ and the phase 
4>i( r / r i) °f the Friedel oscillations in the regime I. 

In terms of variables u and v in the parametric space 
Fig. [21 the condition e r = 1 can be presented as 



1 



where 



r 



(5.8) 



The dependence Eq. (|5.8|) is shown in Fig. 2 with a 
dashed line within the regime I. To the left of this line, we 
have e r < 1, so that 1/r 2 decay of the Friedel oscillations 
is unchanged in the random field. To the right of the 
dashed line, e r is bigger than 1. Then, the dependence 

~Y j 2 

F(r) oc e r , which follows from Eq. (|3.6p . translates 
into faster, but still power-law decay, oc 1/r 7 / 2 , of the 
Friedel oscillations. Note also, that the phase of the os- 
cillations also changes as e r crosses over from small to 
large values. Indeed, as follows from Eq. (|3.6p . we have 
4>(r) -> -7r/4+ l/(2r 3 ) in the limit e r > 1. 
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2. Periodic Magnetic Field 

Consider a particular case of a spatially-periodic mag- 
netic field h{x,y) = h cos(qx). For small enough q the 
"local" description applies. The corresponding condition 
reads 



q<< P° = k *(j£k) 



2/3 



(5.9) 



Under this condition, the averaged Friedel oscillation can 
be found by averaging 

f ~ \ 2 / 3 

Eq. (|3.1[) , in which po is replaced by po I h/ho ) , over 

the distribution, P(h), of the local values of magnetic 
field rather than over the gaussian distribution Eq. (15.5[) . 
This distribution has the form 



P{h) 



(5.10) 



so that instead of Eq. (|5.5[) we have 



oscillates with "period" much larger than the de Broglie 
wave length, but much smaller than the period, 1/q, of 
change of the magnetic field. 

Note that this effect provides a unique possibility to 
measure experimentally the amplitude of a periodic mod- 
ulation. The reason is the following. The envelop 
Eq. (|5.14p due to periodic magnetic field (or electric 
field, i.e., due to the lateral superlattice) translates into 
a distinct low-frequency behavior of the tunnel density 
of states. Namely, the tunnel density of states would ex- 
hibit an "oscillatory" behavior with a "period" u> ~ pov F . 
This period in to depends only on the magnitude of the 
modulation, ho, but not on the spatial period of modu- 
lation, 2tt j q. Therefore, the magnitude of modulation, 
which, unlike the period, is hard to measure otherwise, 
can be inferred from the bias dependence of the tunneling 
conductance. 



B. Friedel oscillations in a random magnetic field: 
Regime II 



— [ dx 

7T J 

-1 



cos(e r x 2 + j3) 



Jo (e r /2) cos ( y + [j 



= Ti(e r ) cos/3 — T 2 (e r ) sin/3, (5.11) 
where J is the Bessel function, e r = {por) 3 /12, and 

(5 ' 12) 

T 2 (e r ) = J ' 



T x (e r ) = J ( y j cos ( y j , 



, sm , 
2 V 2 



so that in a periodic field, instead of Eq. (|3.6p . we have 



2/; 



4>{r) = — arctan 



nl/2 






*(t) 



(5.13) 



t 2 (e r ) 



It is instructive to present the results Eq. (]5 . 13[) in a 
different form, by simply showing how the Friedel oscil- 
lation Eq. p.ip gets modified on average in the presence 
of a periodic magnetic field. Substituting Eq. (|5.13[) into 
Eq. (|3~5)) we get 



Vn{r)) = 



v gV(2k F 



x sm 



2?rr 2 
2k F r - 



Jo 



p3 r 3 

~2A 



(Por) 3 
24 



(5.14) 



Eq. (|5.14[) is a quite remarkable result. It suggests that, 
due to the periodic smooth magnetic field, the averaged 
Friedel oscillations do not get smeared. Rather they ac- 

oscillatory envelope, Jq f%j-Y This envelope 



quire an 



As the magnitude, ho, of the random field decreases, 
the character of semiclassical motion changes from arc- 
like (regime I in Fig. [1} to the snake-like (regime II in 
Fig. HJ) . To estimate for the "widths", Sy, of the snake- 
like trajectories, we use Eq. (|4. 14|) and set x ~ £. This 
yields 



5y_ ^ eh £, ^ ( j_ 
£ mcv F \ k F £ 



1/2 



(5.15) 



Since k F t; 3> 1 and e -C 1 in regime II, we confirm that 

^ C> i- e -i that the snake is "narrow". 

It is clear that at large enough distances, r, the mag- 
nitude, F(r), of the averaged Friedel oscillations falls off 
exponentially with r. The prime question is what is the 
characteristic decay length. As stated in Section IIIII this 
length, r„, is given by Eq. (|3.9|) . Below we derive this 
length qualitatively, and then establish the form of the 
magnitude, F u (r), as well as the phase, 4> u {r), for the 
average Friedel oscillations within the entire domain of r 
by performing the functional averaging of exp(2i<5^o^r)- 



Qualitative consideration 



To recover qualitatively the scale r„ from Eq. (|4. 1 8|) we 
consider the following toy model. Let us divide the inter- 
val (0, r) into small intervals of a fixed length, £ (overall, 
r/£ intervals). Assume now that the random field takes 
only two values, ho and —ho, each with probability, 1/2, 
within a given interval, £. Under this assumption, we find 
from Eq. (|4.13|) A(r) = ho [m(r) — n(r)] £, where m(r) 
and n(r) are the numbers of small intervals within the 
length, r, with h = ho and h = —ho, respectively (obvi- 
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ously, m + n = r/£). From Eq. (|4. 1T|) we get for Sipo- 



hie 



$ fc F | JO 



dx\m(x) — n(x)] ' 



dx\m(x) — n(x)\ 



(5.16) 



Second term in Eq. (|5.16[) is square of the difference (m) — 
(n) of coordinate (not statistical) average values of m(x) 
and n(x). Rewriting m(x) as (m) + Sm(x) and n(x) as 
(n)+Sn(x), and taking into account that 5m(x)+8n(x) — 
0, one can cast Eq. (|5. 16(1 into the form 



Stp, 



•2^2 r r 



O^r 



llx ' dx [5m(x)\ l 



*0 fc F JO 



(5.17) 



Since the typical value of [<5to(x)] 2 is (m(x)) = x/2£, we 
arrive at the following estimate S(po^ r ~ h^r 2 /^k F . 
Equating this additional phase to unity yields r — 
^o^y 2 /hoS, 1 / 2 1 which coincides with r„ defined by 
Eq. (|3.9p within a numerical factor. 



2. Evaluation of the functional integral 

Below we present the analytical derivation of 
Eqs. (|3.7p . (|3.8p . The averaging of exp (2z<5y>o->r) re- 
quired to calculate F u (r), and <^ n (r) from Eqs. ()5. 1|) . (|5.2j) 
reduces to the functional integral 



/ J){fc(r)} exp 


2iSip(r) -W{h(r)} 


fD{/i(r)}exp 


-W{h(r)}~ 





(5.18) 



where 5(p(r) — S(po^ r is given by Eq. (|4. 1T[) . and 
exp(—W{h}) with W {h(r)} given by 



W{h} 



dx\ dyi 



dx 2 dy 2 



5 ''0 JO J -co JO 

x h(xi,yi)h(x 2 ,y2) - £2,2/1 - J/2), (5.19) 



is the statistical weight of the realization, h(x,y). The 
dimensionless function «(r, r') is related to the correlator 
Eq. (|2.ip in a standard way 



/" dr'K(r, r')K(r', r") = £ 4 <5(r - r"). 



(5.20) 



The reason why the functional integral Eq. (|5.18|) can 
be evaluated explicitly is that both W{/i} and 5<po^r 
are quadratic in the random field, h(x,y). The fact that 
we integrate over realizations of h(x, y) defined on the 
interval which is finite, < x < r, in the x-direction and 
infinite in the y direction suggests the following expansion 
of h(x,y) 



00 pQO 

h(x, y) = h ^ / dqA n ,q( 

J -co 



exp 



^2irinx\^ 



(5.21) 



The asymmetry between x and y is quite significant in 
the calculation below, namely, for r ^> £, the charac- 
teristic values of x turn out to be much larger than the 
characteristic values of y is ~ Sy <C £, see Eq. (|5.15l) . 
This allows to replace K(a;, y, x',y') in Eq. (|5.20p by 
7£ K(0, y — y')S(x — x'), where the dimensionless con- 
stant 7 is defined by the relation 



00 00 



J dx J dy K(x,y) 
7 = ° Z2E = (1) 

I 00 In/ 

C/dyK(0,y) 



J dz zK-(z) 


00 

/ dz K ( z ) 




-, (5.22) 



where in the second identity we used the fact that K(x, y) 
is isotropic. Substituting Eq. (|5.2ip into Eq. (|5.19[) . we 
obtain 



00 

w w = ^ E /* 



I , q I 



(5.23) 



where K.(q) is the Fourier transform of the correlator, 
more precisely, 



K(q) 



dy 



iQv/6 



2ir J £ 



K(0,y). 



(5.24) 



Expression for 5ip(r) in terms of the coefficients, A n ,q, 
follows upon substitution of Eq. (|B3|) into Eq. (|4.18|) 



= EE/ d 9lAu,gi / dq 2 A n2 , q . 

F ni=— 00 n 2 = -oo 

/ dxi dx 2 \r ^—^- — x 2 Q(x 2 — Xi) 

J0 Jo 1 T 

/ 2?ri 



exp < (ni.Ti + n 2 X2) 

I r 



—x\Q(xx - x 2 ) 
Performing the integration, we obtain 



(5.25) 



MO = 5" I / *Z ^0,, + E C »| / d 1 A 
^ { J n>0 J 

+ Jdq A , q J dqJ2[b 



(5.26) 



where numerical coefficients b n and c n are defined as 

1 



1 i 
+ 



2n 2 n 2 2ixn 



2it 2 n 2 ' 



(5.27) 



In writing the result of integration in the form Eq. (|5.26D 
we have used the dimensionless parameter e defined by 
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Eq. (|3.10p . The meaning of this parameter is the addi- 
tional phase Eq. (|4.17[) . acquired by the electron travel- 
ling the distance ~ £ in a constant magnetic field, ho. 
Since our calculation pertains to the limit r 3> £, the 
relevant values of e are small. 

The functional integration reduces now to the infinite 
product of the ratios of integrals over A n , q and A* x q . The 
details of calculation are given in Appendix [Cj Here we 
present only the final result for r>( 



1 



(e 2td ^ r >) = - 

i 2i ( r 



3 V r n 



n 



L n 2 — 2i(r/r n ) 2 /TT 2 ' 
where the characteristic length, r„, is defined as 



2i 



4fc F $ 2 1 



Va y ( 27r )i/2 7 ^2- 



(5.28) 



(5.29) 



The above definition specifies the numerical coefficient, 
7], in Eq. J331) of Section HH as tj = 2/(2vr) 1 / 4 7 1/2 - This 
coefficient depends on the explicit form of the correla- 
tor via the factor 7, given by Eq. (|5.22[) . It is seen that 
r„ ~ ^/e 1 / 2 is indeed much larger than £. This means 
that, in the regime II, Friedel oscillations survive well 
beyond the correlation radius of random magnetic field. 
Note also a distinctive dependence r u oc 1/ho of the char- 
acteristic scale on the magnitude of the random field. In 
fact, the infinite product in Eq. (|5.28|) can be evaluated 
for arbitrary r/r n , using the identity 



SIM 



which yields 



2iS<p(r)\ _ 



1 



(l + i)(r/r n ) 



(5.30) 



(5.31) 



sin(r/r„) cosh(r/r n ) + i cos(r/r„) sinh(r/r n ) 



With the help of Eq. (|5.31[) we can calculate the magni- 
tude, F n (r), and the phase, (/> n (r), of the Friedel oscil- 
lations in the regime II. Corresponding expressions are 
given by Eqs. fl37|) and ([378}) . 



3. Limiting cases 

It is not surprising that Friedel oscillations in the 
regime II are smeared more efficiently than in the regime 
I. The small-r and the large-r asymptotes of F u (r) are 
the following 



45 V rn 



r <C r ls 



(5.32) 
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FIG. 6: (Color online) Third-order process describing creation 
of a pair by initial electron at point r = 0, rescattering within 
the pair at point r = ri , and annihilation of the pair at point 
r = Y2- Diagram corresponding to this process is shown later 
in the text (first diagram in Fig. [8]). 



Fn(r) =3V2pi)exp(- 



T 



r>r„. (5.33) 



We see from Eq. (|5 . 33[) that Friedel oscillations decay ex- 
ponentially as r exceeds r n . This should be contrasted to 
Eq. (|3.6p for the regime I, where the Fj (r) falls off slowly, 
as r~ 3 / 2 , with r. On the qualitative level, the strong dif- 
ference between the regimes I and II, that is reflected in 
the different characters of decay of and F u (r), is 

that in regime I the random field does not change within 
the characteristic spatial interval, r l7 while in regime II 
the sign of the random field changes many times within 
the characteristic spatial interval, r„. 



VI. DENSITY OF STATES: QUALITATIVE 
DISCUSSION 

In the previous consideration we had demonstrated 
that in two regimes of electron motion in random mag- 
netic field, i.e., regime of arcs, I, and regime of snakes, 
II, there are two length-scales, Ti and r n , respectively 
that govern the interaction effects. In this section we 
demonstrate that the density of states, 8v(oj), exhibits 
an anomalous behavior within the frequency range u> ~ 
v F /r 1 in the regime of arcs, and oj ~ v F /r n in the regime 
of snakes. 

The process underlying the interaction corrections to 
the density of states is creation (and annihilation) of the 
virtual electron-hole pairs by an electron moving in the 
random field. Our central finding is that, unlike the case 
of point-like impurities^, the low-w structure in the den- 
sity of states emerges as a result of electron-electron scat- 
tering processes involving more than one pair. 

We start with a three-scattering process in the regime 
of arcs, and demonstrate qualitatively how the frequency 



scale, v w /rj, emerges. Three-scattering process involves 
two virtual pairs. Consider first this process in the ab- 
sence of the random field. It is illustrated in Fig. [6] In 
analysis of this proces a 36 i 65 ' 66 it was established that the 
directions of momenta of the participating electrons are 
strongly correlated, namely, they are either almost par- 
allel or almost antiparallel. Quantitative estimate for the 
degree of alignment of the momenta can be obtained from 
inspection of Fig. [6] If the scattering acts take place at 
points 0, r l5 and r 2 , then the corresponding matrix ele- 
ment contains a phase factor 



exp [2ik F (n - r 2 + \r 1 - r 2 |)] . (6.1) 
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FIG. 7: Diagrams for the second-order corrections Eq. (|7.2p 
(left) and Eq. (|7,3[1 (right) to the density of states. 



VII. DENSITY OF STATES: ANALYTICAL 
DERIVATION 

A. Absence of a zero-bias anomaly in the second 
order in the interaction strength 



This phase factor does not oscillate, if the angle between 

1 /2 

the vectors ri and r 2 is smaller than (l/fc F r) ' , where r 
is the typical length of , r 2 . 

The above angular restriction constitutes the origin of 
a zero-bias anomaly in the regime of arcs. Zero-bias 
anomaly emerges as a result of the suppression of the 
three-scattering process in the field, h . This suppres- 
sion is due to curving of the electron trajectory by the 
angle ~ r/R^, see Fig. [5j and it occurs when the curving 
angle exceeds the allowed angle of alignment. Therefore, 

1/2 

upon equating (l/fc F r) to r/i? L , we find r = r 11 which 
leads us to the conclusion that u> ~ fp/fi is the energy 
scale at which 5v{ijj) exhibits a feature. Note that, in 
considering the Friedel oscillations, we inferred the scale 
Ti from a different condition, namely, that the additional 
phase, ~ (po r i) 3 i due to the elongation of a trajectory 
in magnetic field is < 1. Thus we conclude that, in the 
regime of arcs, the same spatial scale, r u which governs 
the "dephasing" of n 2 fc F (r) (a polarization bubble) also 
governs the suppression of the three-scattering process, 
which involves three loops . 

The above analysis of phases in the matrix element 
of the three-scattering process can be extended to the 
regime of snakes. This analysis yields that three- 
scattering process is efficient at distances r < r„, see 
Eq. (|3.9p . Analysis of phases similar to the phase, given 
by Eq. (|6.1[) . also suggests that two-scattering processes 
are insensitive to the magnetic field. This insensitivity 
can be explained as follows. Calculation of the contri- 
bution to the density of states from the three-scattering 
process with matrix element Eq. (|6.ip involves integra- 
tion over positions of ri and r 2 , with respect to the ori- 
gin, r = 0, which reveals the angular restriction on their 
orientations. Similar integration for a two-scattering pro- 
cess involves only the orientation of the interaction point, 
r, with respect to the origin. Then the angular restric- 
tion, and its lifting by magnetic field, does not emerge. 
In the next subsection the above qualitative arguments 
are supported by a rigorous calculation. 



We start from general expression for the average den- 
sity of states 

5v{w) =-- (im G a) (r,r)\ , (7.1) 

7T \ / h(x,y) 

where (...) denotes disorder averaging defined by 
Eq. (|5.18p . In the second order in interaction strength, 
the random-field-induced correction to the density of 
states are determined by two diagrams shown in Fig. [7] 
The corresponding analytical expressions read 

2 f dSl f 

5vi(u) = 41m- 2~ dr dr 1 dr 2 G u (r,rt) 

x Go(r 1 ,r 2 ){l/ 2 (2fc F )n 2feF (r 1 ,r 2 ,c i ;-r!) 
+ F 2 (Q)n (r 1 ,r 2) c i ;-O)}G w (r 2l r) I (7.2) 

2 f dfl f 

6v2(u) = -2 Im- / —— I dr dr 1 dr 2 G u (t,ti) 

x G n (ri,ra)G u (r 2 ,r){7(0)[2V(2A! P ) - V(0)] 

x n 2feF ( ri , r 2 , u - O) + V 2 (0)n (ra , r 2 , u> - fi)} , 

(7.3) 

where V(0) and V (2fc F ) are the Fourier components of the 
interaction potential V(r) with momenta zero and 2fc F , 
respectively. Three Green functions in Eqs. (|7.2[) . (|7.3|) 
describe the propagation of electron between the points 
(r,ri), (ri,r 2 ), and (r 2 ,r), Fig. [7] . Polarization bub- 
ble describes the creation of electron-hole pair at point 
ri and annihilation at point r 2 . Difference in signs in 
Eqs. (|7.2|) . (|7.3p is due to the fact that the first diagram 
contains two closed fermionic loops, whereas the second 
diagram contains only one. Numerical factors 4 and 2 
in Eqs. (|7.2[) . (|T.3[) come from summation over the spin 
indices. The difference between them is due two the fact 
the spin of electron-hole pair is not fixed in the first di- 
agram, but it is fixed in the second diagram. The factor 
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FIG. 8: Third-order diagrams contributing to the zero-bias 
anomaly in the density of states. Random field enters via the 
phases of the Green functions. 



2 in the product 2V(0)V{2k F ) in Eq. (JTrSJ) is related to 
the annihilation of the electron-hole pair, since the hole 
is annihilated with initial electron. Then the momentum 
transfer can be 2k F in the course of creation and zero in 
the course of annihilation, and vice versa. 

It is important to emphasize that the Green functions 
and polarization operators in Eqs. (|7.2p . (|7.3p contain the 
information about the random field, h(x,y), via their 
additional phases: (^ ri ^ r2 in G u {ri,r 2 ) and 2ip ri ^ r2 
in Il2fc F (ri,r2, w). The phase, ip ri -, r2 , always enters in 
combination with a main term, fc F |ri — r2 1 . Obviously, 
IIo (i"i , r2 , o->) does not contain a field- induced phase. 
Thus, only the terms containing U 2 k F in Eqs. (|7.2[) . (|7.3p 
should be considered. 

Now it is easy to see that Svi and 8v2 do not ex- 
hibit a field- induced anomaly at small to. This is because 
the field dependence is cancelled out in the integrands of 
Eqs. (|7.2[) . (|7.3[) . To see this, we first note that the inte- 
gration over r in Eqs. (|7.2p . (|7.3p can be easily performed 
using the fact that J drG^(ri, r)G u; (r, r 2 ) is equal to the 
derivative, 3G UJ (ri,r2)/duj. Then we note that the con- 
tribution to 8v\, 5v2 comes only from "slow" terms, in the 
product of two Green functions, G w (ri,r 2 ), Gn(i"i,r 2 ), 
and Il2fc F . These slow terms do not contain rapidly os- 
cillating factors exp{2ifc F |ri — r 2 |}. On the other hand, 
cancellation of the rapid terms in the product automat- 
ically results in the cancellation of the field-dependent 
terms. 

As it was explained in qualitative discussion, the situ- 
ation changes in the third order in the interactions. Cor- 
responding expression for 8v{lj) is derived in the next 
Subsection. 



B. General expression for the third-order 
interaction correction to the density of states. 

Relevant diagrams for the third-order correction to the 
density of states are shown in Fig. [5] The same 8 dia- 



grams were considered in Ref . |36j in the momentum space. 
In Ref. l36| the analysis of these diagrams was restricted 
to small momenta. In our coordinate representation this 
means that only Ilo(r) parts of the polarization opera- 
tors was kept, whereas Il2fc F (r) parts were neglected. As 
explained above, to reveal the sensitivity to the random 
field, we will keep only the Il2fc F (r) parts. Then the cor- 
rection to the Green function corresponding to the sum 
of eight diagrams in Fig. [8] acquires the form 

5v{uo) = 2V(0)V(2k F ) \2V{2k F ) - V(0) (7.4) 

x Im ^2" J ^~ J dr dr i dr 2 G w (r,rx) G (ri,r 2 ) 
x n 2 fc F (ri,0, oj — ft) n 2 /c F (0,r 2 , a; — fl) G 

All the diagrams reduce to the same integrals. Concern- 
ing the difference in numerical coefficients, it comes from 
the number of closed fermionic loops and the spin de- 
grees of freedom. Taking this into account interaction 
coefficient corresponding to the first two diagrams will 
be 2 • (-2) 2 V 3 (2k F ). Coefficient of the third diagram 
is (— 2) 3 V 3 (2k F ). Thus we see, that the contributions 
tx V 3 (2k F ) cancel each other. 

The first and the second diagrams in the second row 
are equal to each other, and each of them has a coefficient 
(~2) 2 V(0)V 2 (2k F ). Coefficient of the last diagram in the 
second row is (— 2)V 2 (0)V (2k F ) , since it has only one 
closed fermionic loop. Finally, the first diagram in third 
row has only one closed fermionic loop and is equal to the 
second diagram on the third row. Each of these diagrams 
contributes with the coefficient (— 2)^(0) V 2 (2fc F ). 

On the physical level, 8 diagrams in Fig.[5]describe dif- 
ferent electron-electron three-scattering processes. For 
example, the first diagram corresponds to creation of 
electron-hole pair by the initial electron followed by 
rescattering within a created pair and, finally, its anni- 
hilation. Three stages of this process are illustrated in 
Fig. [5J However, creation, rescattering, and annihila- 
tion of a pair can follow a different scenario, namely, the 
rescattering process can involve the initial electron. This 
scenario is captured by the second diagram in the first 
row in Fig. [H] 

At this point, we note that diagrams in Fig. [5] do not 
exhaust all possible three-scattering processes. In fact, 
all diagrams in Fig. [8] have identical structure, in the 
sense, that they can be combined into a single generalized 
diagram, as shown in Fig. There are also eight other 
diagrams combined into a single generalized diagram, as 
shown in Fig. [9Jd that are not sensitive to the random 
field. This is because, in the absence of the random field, 
the phase factor corresponding to Fig.[3jD is large, namely, 
2-2fc F (n +r 2 ). 

The crucial difference between the contributions 
Eqs. (j7T2|) , (f773|) and Eq. (fT4|) is that the cancellation of 
the rapid-oscillating terms in in the integrand of Eq. (|7.4p 
preserves the field-dependence. To see this, we first re- 
place / drGo(ri, r)Go(r, r 2 ) by <9Gn(ri, r 2 )/dtt, as dis- 
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FIG. 9: (a) Eight diagrams for <5G^(r, r), that are shown in 
Fig. [8] are combined into one generalized diagram. Electron- 
electron scattering processes take place at points 0, ri, and r2\ 

(b) Eight third-order diagrams that do not contribute to the 
zero-bias anomaly are combined into one generalized diagram; 

(c) Two types of four-leg interaction vertices are combined 
into big dots; (d) An example of a third-order diagram of 
type (b). 



cussed above, and then consider the phase of the product 

Gn(ri,r 2 ) G^(r 1 ,r 2 ) 

xlW0,r 2 ,w-r!) n 2 fc p (ri > 0,u;-n). (7.5) 

Fig. [SJs illustrates this product graphically. It is seen 
from Fig. 03 that, when the fast oscillating terms 
exp{2i/c F |ri — r2 1 }, exp{2ifc F ri}, and exp{2i/c F r 2 } can- 
cel each other out, the additional phase enters into the 
product either in combination 



or in combination (see Fig. 03 ) 



(7.6) 



(7.7) 



Since additional phases defined by Eqs. (|4.18[) . (|4.19|) 
are cubic in distance, the combinations Eq. (|7.6p and 
Eq. (|7.7p are nonzero. This is in contrast to the two- 
scattering processes, where the cancellation occurs iden- 
tically for arbitrary dependence of 6ip(r) on r. In turn, 
non-cancellation of additional phases in Eqs. (|7.6p , (|7.7p 
means that the random field causes a zero-bias anomaly, 
more specifically, a feature in 5v(u>) at small w. 

The final form of 6v(u)) emerges upon integration of 
Eq. (|7.4p over azimuthal angles of ri and r 2 , which can 
be performed analytically, using the relation 



ip(ri+r 2 ; 



sin [p (ri ± r 2 ) + 7r/4] 



(7- 



Upon combining rapidly oscillating terms in the inte- 
grand of Eq. (|7.4[) into "slow" terms, we obtain 



Sv(u) = 5v {+) {lo) + 8v { -\uj), 



(7.9) 



where 



dr\dr 2 

W 2E ¥ ^/ 2 k 1 J 2 Jr 2 >r 1 {rir 2 f /2 

x^+ra) 1 / 2 / dnsm[v- 1 (ij-n)(r 1 +r 2 )] 
Jo 

x sin +1- t±3 {ri + r A , (7.10) 

4 -y F 



and 



8y{-\w) _ (y Vf f dndr 2 
2E F ir 3 / 2 kl /2 Jr 2 > ri {rir 2 ) 3 / 2 

x(r 2 -ri) x / 2 / dfl sin [v- 1 (w-fl)(r 1 +r 2 )] 



x sin 1 2*4 ] + \ + ^7^( r2 " ri )} ' ( 7 ' n ) 

where we had assumed that the interaction is short- 
ranged and set V(0) = V(2k F ) = i>oV ■ Two contri- 
butions in Eq. (|T.9[) correspond to the locations of the 
points ri and r 2 on the opposite and the same sides from 
the origin, respectively, see Fig. [SJs. 

We note that the phases 6ip^~\ 8<p^~\ which enter into 
the argument of sine in Eqs. (|7.10[) . (|7.1ip . are quadratic 
in the random field, h(x,y), as seen from Eqs. (|4.18D . 
(|4.19[) . This suggests that the averaging over realiza- 
tions of h(x, y) can be carried out analytically in the 
integrands of Eqs. (|7.10|) . (17. lip . Similarly to the case 
of Friedel oscillations, it is convenient to perform this 
averaging separately for the regimes I and II. This is 
done in Sections IVIIIi IIXI below. In the remainder of 
this Section we will evaluate the interaction correction, 
5v{u>), for two particular cases: (i) constant magnetic 
field, h(x,y) = ho, in a clean electron gas, and (ii) 
h(x, y) = ho in electron gas with small concentration of 
point-like impurities. 



C. Case of Constant Magnetic Field: Oscillations 

of Su(co) 

In a constant magnetic field h(x,y) = ho the char- 
acteristic scale of frequency in Eqs. (|7.10[) . (|7.11[) is 
luo = v F /r l . This was stated in Section IIIH Now this 
scale of frequencies emerges naturally upon substituting 



in Eqs. (f77TO|) . J77TT]) the phases 25<p\ T> , 28<p\ T> , calcu- 
lated from Eq. (14. 1 8[) in a constant magnetic field 



3 

28ip < £ ) =T^r l r 2 (r 1 ±r 2 ) 



(7.12) 

where po is defined by Eq. (|3.2p . The integrals in 
Eqs. (|7.10p . (|7.1ip converge at distances r\,r 2 ~ p^ 1 = 



15 



2 

Si 



o 
-1 



1 1 1 


1 


1 


- 


\ 








■ 


\ 










\ 








- 










■ 


















X 




i 


■ 


■ 


■ 





I i I i I i L 

12 3 



FIG. 10: (Color online) Dimensionless correction Eq. (|7. 13[) 
to the tunnel density of states in a weak constant magnetic 
field is plotted vs. dimensionless energy x — 2 2 / 3 lu/ujq. The 
plot is obtained upon numerical integration in Eqs. (|7.10|l . 

CllJJ. 



prefactor in Eq. (|T. 13|) that the magnitude of 5v scales 

3/2 

with ho as r r oc ho- We emphasize that the correc- 
tion 8v(w) remains distinguishable even when the struc- 
ture in the density of states due to the Landau quanti- 
zation is completely smeared out, e.g., due to finite tem- 
perature. This follows from the above relation between 
loq and the cyclotron frequency, uj c , namely, (u c /luq) ~ 
(uj c /E f ) 1/3 « 1. 

In discussing the relevance to the experiment one 
should have in mind that realistic samples always con- 
tain certain degree of disorder. Therefore, the question 
remains as to whether the oscillations of Su(u>) in a con- 
stant magnetic field survive in the presence of the short- 
range impurities. This question is non-trivial, since im- 
purities themselves give rise to the singular correction 
to Sv(u>) (zero-bias anomaly) even in a zero field. Then 
the above question can be reformulated as: whether the 
field-induced oscillations are distinguishable on the back- 
ground of a zero-bias anomaly. It turns out that, by 
introducing the Friedel oscillations, point-like impurities 
actually enhance the oscillatory part of 5v{uj). This ques- 
tion is addressed in the next subsection. 



tv As a result, 5v^ + ' and are certain universal 

functions of ujr 1 /v- e — lo/ujo- The plot of + 5v^' 

vs. dimensionless ratio x = 2 2 ^lj/ujq is presented in 
Fig. [TDJ To isolate the frequency dependence, in addi- 
tion to x, we had introduced the dimensionless variables 
r-i/r-i and r 2 /r 1 after which 5v{uj) acquires the form 



Su(oj) 

V<3 



2 2 / 3 ( 7 rfc F r I ) 3 / 2 



B(x). 



(7.13) 



The integral over f2 in Eqs. (|7.10|) . (|7. 1 1[) can be eval- 
uated analytically. The remaining dimensionless double 
integrals were calculated numerically. While the char- 
acteristic scale, x ~ 1, of change of the function B(x) 
follows from qualitative consideration, Fig. [TO] indicates 
that B(x) also exhibits sizable oscillations. These os- 
cillations come only from the contribution . They 
owe their existence to the peculiar structure of the ar- 
gument of sine in Eq. (|7.10p . Namely, this argument 
has saddle points with respect to both r\ and r 2 at 
n = r 2 = 2 1 / 3 r I (ti>/aj ) 1 / 2 /3 1 / 2 Oscillatory behavior of 
B(x) is governed by the value of the argument at the 
saddle point, which is ~ (uj/ujo) 3 ^ 2 - Strictly speaking, 
the saddle point determines the value of the integral only 
when uj 3> loq. However, numerics shows that oscillations 
in Fig. [TOJ set in starting already from x ~ 1. These os- 
cillations reflect the distinguished contribution from the 
three-scattering process, shown in Fig.03, in which scat- 
tering events occur at n = r 2 = 2 1 / 3 r I (w/o;o) 1 ^ 2 /3 1 / 2 . 

Eq. (|7.13[) and Fig. [TO] constitute an experimentally 
verifiable prediction. Correction Eq. (|7. 13[) describes 
the the feature in the tunneling conductance of a clean 
two-dimensional electron g clS £1S 0j function of bias that 
emerges in a weak magnetic field, ho- It follows from 



D. Ballistic Zero-Bias Anomaly in a Constant 
Magnetic Field 

Conventional ballistic zero-bias anomaly 3 -, caused by 
point-like impurities, is described by two second-order 
diagrams, shown Fig. in which one of two inter- 
action lines is replaced by an impurity line. As was 
shown in Ref. U these diagrams with one interaction 
line and one impurity line yield a singular correction, 
8u{uj)/vq ~ (v V/E f t) ln(w), to the density of states. 
Here r _1 is the scattering rate proportional to the impu- 
rity concentration. Qualitatively, the singular correction 
originates from the combined scattering of electron by 
the impurity and the Friedel oscillation oc sm(2fc F r)/r 2 , 
created by the same impurity. This Friedel oscillation 
is represented by the polarization loop in Fig. [7] In 
the presence of the impurity, this loop describes static 
response of the electron gas, and thus the polarization 
operator, Tl 2 k F (uj , r) , corresponding to the loop should 
be taken at u> = 0. As was mentioned in Section IIII1 
a weak perpendicular magnetic field, h, leaves the loga- 
rithmic correction unchanged. To reveal the sensitivity 
to h, one should calculate 8v to the next (second) or- 
der in V . Corresponding diagrams with one impurity 
and two interaction lines are shown in Figs. 111! 1121 and 
ITOl It is easy to see that there are overall 24 different 
diagrams. Indeed, the generalized diagram, Fig. [9ja), 
for the third-order interaction correction contains three 
generalized four- leg vertices shown in Fig. [He). Hence, 
Fig. [5] (a) represents 2 3 = 8 different diagrams. In each 
of these 8 diagrams, the impurity line can replace inter- 
action line in three places, generating one of 24 different 
diagrams that are shown in Figs. [Til [TOJ andfTOl All these 
diagrams are divided into three groups according to their 
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FIG. 12: 8 out of total 24 second-order diagrams for ballistic 
zero-bias anomaly in a weak constant magnetic field, which 
contain one dynamic polarization operator. 



1 



FIG. 11: Second-order diagrams contributing to the oscillat- 
ing part (see Fig. [4| of the ballistic zero-bias anomaly in a 
weak constant magnetic field. Magnetic field enters through 
the phases Eq. (|4.7|l of the Green functions. Dashed line repre- 
sents the impurity scattering. All 12 diagrams (a)-(l) contain 
two static polarization operators. 



dependence on tj. Namely, all 12 diagrams in Fig. fTTl havc 
the same ^-dependence. Similarly, the w-dependence of 
all 8 diagrams in Fig. [T2] is the same. This also applies 
to 4 diagrams in Fig. 1131 However, the corresponding lu- 
dependencies are slightly different from each other. The 
origin of this difference can be traced from comparison 
of diagrams Fig. QT] (a), Fig. HU (a), and Fig. [TJ (b). 
Diagram Fig. [IT] (a) contains two polarization loops sep- 
arated by the impurity line. As a result, the expression 
corresponding to this diagram, contains two static polar- 
ization operators, n 2 fc F (0, r). Diagram Fig. [T3] (a) con- 
tains one finite-w polarization loop, H 2 k F {<-^ , r) . Finally, 
the diagram Fig. [T3] (b) does not contain polarization 
operators at all, but rather a different object, namely, a 
polarization loop crossed by the impurity line. Important 



is that the expression, corresponding to this object 

JJ( w -n,| ri -r 2 |) = -i J ^L Gai (0,'i) 

x G w -n+fii (ri , 0)G a ,_ n +n 1 (0, r 2 )G 0l (r 2 , 0), (7.14) 

contains a "fast" part, ri2fe F ( CJ ' 7 ')' which oscillates as 
exp(2i/c F |r! — r 2 |), i.e., in the same way as polarization 
operator. 

The full analytical expression corresponding to the di- 
agram Fig. [11] (a) reads 



Svi(u>, h) = Im 



W 2 (2k F ) 



drdridr 2 G w (r,ri) 



xG w (ri,r 2 )n2fc F (0,ri)n2fe F (0,r 2 )G w (r2,r) 
6V 2 (2k ] 



Im- 



J drx dr 2 d u G UJ {Y l ,T 2 ) 
xG w (r 1 ,r2)n2 fcF (0,r 1 )n 2fcF (0,r 2 ), (7.15) 



TT 2 VqT 



where in the second identity we had performed integra- 
tion over r. 

Analytical expression for the diagram Fig. (a) has 
the form 



2V 2 {2k F 



5v 2 (u,h) = -Im , 

2tt z vqt 

x G u (r,n)G u (r u 0)G 



dr dridr2 



(7.16) 



— G o (0, ri)G n (ri,r 2 )n 2fcp (w - SI, \n - r 2 |). 
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FIG. 13: 4 out of total 24 second-order diagrams for the ballis- 
tic zero-bias anomaly in a weak constant magnetic field, which 
contain a polarization loop crossed by the impurity line. 



Finally, the expression for the diagram [T3] (b) is the fol- 
lowing 



8u^{uj, h) 



-Im 



2V 2 {2k F ) 
2tx 2 vqt 



■Jdr dY 1 dr 2 G w (r, ri)G w (r 2 , r) 

* J to J ^7 Gn ( r i' r2 ) Gfi i(°' ri ) 

xG w -n+n 1 (r x , 0)G w _ n+Ol (0, r 2 )G 0l (r 2 , 0). (7.17) 

Upon integration over r, it can be expressed through 
n 2feF M, defined by Eq. (EH, as 

6va(u, h) = — — [dr 1 dr 2 9 u ImG u (ri, r 2 ) 

x/ — ImG n (ri,r 2 )Im TT (w - fi, \r x - r 2 |). 

(7.18) 

Despite all 12 diagrams in Fig. [TT1 have the same fre- 
quency dependence, their prefactors represent different 
combinations of V 2 (2k F ), V 2 (0), and V(2k F )V(0). The 
same applies to 8 diagrams in Fig. 1121 and to 4 diagrams 
in Fig. [T3J Taking into account the numerical factors 
in these combinations amounts to the following replace- 
ments: in 5v\ 



W 2 (2k F ) -> 3F 2 (0), 



(7.19) 



in 5v 2 
-2V 2 (2k F ) 

and in Si's 
-2V 2 (2k F ) 



V(0)V(2k F ) - V 2 {2k F ) - V 2 (0) 



(7.20) 



V(0)V(2k F ) - V 2 (2k F ) - V 2 (0) 



(7.21) 



These replacements must be taken into account when 
calculating the full correction 6v(u>) from 5v\, 8is 2 , and 
Su 3 . 



Below we demonstrate that all three contributions 6vi , 
8v2, and are oscillatory functions of lo. Detailed 
derivation will be presented only for 8v\. 

Analogously to the derivation of Eqs. fTITD]) . (fTTTTjl . 
we can perform the integration over the azimuthal 
angles of ri and r 2 analytically using Eq. (|7.8[) . 
Then, extracting a "slow" term from the product 
of trigonometrical functions, we obtain 5v\{lj) = 

[6^V 2 {2k F )/E F r] (w /^f) 1/2 Pi(2 2/3 w/w ), with lo = 

2 x I 3 lo 2 J 3 E F ^ 3 , where the function Pi(x) is denned as 



where 



p+ f d Pl dp 2 J 



(7.22) 



(Pi+p 2 ) 1/2 



x(pi +Pi)-j- piPi{,Pi + 92) 



cos 



x(pi +p 2 ) - ~ } j 



(7.23) 



pre 



dp\dp 2 



(P2-Pl) 1/2 



p 2 > Pl {PlP2) 3/2 { 
| COS x(p 2 - Pi) + ^ + P\Pl{P2 - Pi) 



(7.24) 



Here the constant factor, tr, is given by er = (S^ 1 / 6 )/^ 3 / 2 . 
In Appendix [El we demonstrate how the function Pi (a;) 
can be cast in the form that is convenient for numerical 
evaluation and extracting asymptotes. This form is given 
by the following double integral 

„ , s f°° dz f° dv , « 

Px x )=Aa / -— / -== 7.25 
Jo z A/ J _ 4 + 4 







7T 


z 3 ' 






7T~ 




cos 


ccz - 


"4" 




V 


— COS 


xz - 


" 4- 


) 



The fact that Pi(x) oscillates at large x 3> 1 follows 
from the observations that (i) first cosine in the brackets 
in Eq. (j7T25)) has a saddle point z = (x|i;|/3) 1/2 , and (ii) 
the major contribution to the integral over v comes from 
the lower limit v = —4 (corresponding steps are outlined 
in Appendix [Dj) . This yields 



Pi(x) 



25/339/4 
TT 1 / 2 



1 



(I) 



3/2 



(7.26) 



The argument x 3 / 2 in the cosine in Eq. (|7.26j) can be 
presented as w 3 / 2 / {2}/ 2 Ljj c El^ 2 ^j , so that the "period" 
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in to is much bigger than the cyclotron energy, u) e , as was 
discussed above. 

The analysis of the contributions bv 2 (w) and bv^iu) 
can be carried out in a similar way. They exhibit the 
same oscillations as Eq. (|7.26p . The difference is that, 
due to integration over O in Eqs. (|7.16[) and (|7.18[) , both 
5v 2 (w) and 5vs(lo) contain an extra factor lo/loq, see 
Eq. (|3. 13|) . and thus their contribution to the net cor- 
rection 8v is dominant at lo ujq. 



where the functions S+, <SL, C+, and C_ are defined as 



S±(pi,p 2 ) = {pi±p 2 ) 1/2 sm 7T(z+ z')(fii ± P2) 



x i Ui[ p 1 p 2 (pi ±p2j ) -Vn>, 



(8.7) 



C±( Pl ,p 2 ) = (pi ±p 2 ) 1/2 cos '- T {z + z'){pi±p 2 ) 



*U 2 p\p 2 {p\ T p 2 ) 



VIII. ZERO-BIAS ANOMALY IN THE 
AVERAGED DENSITY OF STATES IN REGIME I 



With the help of the identity Eq. (|5.5p the integrand 
in the average Sv(u>) can be expressed in terms of func- 
tions Ui_ 2 [rir2(r 2 ±ri)po/4], where the functions U\_ 2 
are defined as 



Ui{x) 




U 2 (x) 




1 + x 2 : 



(l+x 2 ) 1 / 2 -! 
1 + x 2 



(8.2) 



Upon introducing dimensionless variables p\ 2 
Po r i, 2 /2 2 ^ 3 , we present the final result in the form 



(8.3) 



with 



^0 = v F p = 2E F 



ho 



and with constant, C, defined as 
(yoV) 3 ( i>. 



2/3 



(8.4) 



C 



2tt Vfcp$o 



»''?)*'. (8.5) 
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The dimensionless function, T{z), describing the shape 
of the anomaly, is given by the following double integral 
over pi , p 2 



l(z)=2+(z)+l-(z) 



(8.6) 



dpidp 2 
(Pip 2 ) 3/2 

P2>pl 



dz sin 



(z- z')(pi +p 2 ) 



[s + (pi, p 2 )+C + (pi, p 2 ) + S_(pi, p 2 )+C-(p 1 , p 2 )} , 



In definitions of S + and S- we had subtracted from the 
function Ui(a) the zero-field value Ui(0) = ^/tt. Inte- 
gration over z' in Eq. (|8.6p can be easily carried out 
analytically. The remaining integrals over p±, p 2 were 
evaluated numerically. Direct numerical integration en- 
counters difficulties due to very fast oscillations of the 
integrand in Eq. (|8.6[) . These difficulties can be over- 
come by a proper change of variables in the integrand. 
This procedure is described in Appendix D. The result- 
ing shape of the zero-bias anomaly is shown in Fig. [14] 
The small- z <C 1 behavior of T{z) is 81nz, i.e., it diverges 
logarithmically. The cutoff is chosen from the condition 
that T{z) approaches zero at large z. Note, that 1{z) 
exhibits a pronounced feature around z = X. The origin 
of this feature lies in strong oscillations of the integrand 
in Eq. ()7.9|) . The "trace" of these oscillations survives 
after averaging over the magnitude of the random field. 
In fact, the oscillations persist beyond z = 3. This is 
reflected in the z» 1 asymptote of the function I(z), 



l+(z) 



_ sin (2 8 / 3 V3z) 

- 23/4 ^ ^ ,3/4 ' ^ 



{-2 8 / 3 z} 



(8.9) 



To derive this asymptote, it is more convenient to first 
take the limit of large u> in Eq. (|7.9[) and perform the 
averaging over the random field only as a last step. In 
the limit uj » uj following simplifications of Eq. (|7.9|) 
become possible. Firstly, the second term in the square 
brackets can be neglected, since it does not produce oscil- 
latory contribution to bv. Secondly, one can set SI = in 
the integrand, so that the integration over O reduces to 
multiplying by lo. Lastly, upon converting the product of 
sines into the difference of cosines, one finds that the in- 
dependence is present only in the term, corresponding to 
the difference of arguments. As a result, the oscillatory 
part of 8v(u>) at u> ^> loq acquires the form 




(is V) 3 ojujI /2 r dpidp 2 

215/6^3/2^3/2 J ( pi p 2 f/2 
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(Pl+P 2 ) 1/2 



P\Pl{P\ + Pi) + 



2 5/3 



CJ 
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x(Pi + P 2 ) 



(8.10) 
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This cancellation, as we demonstrate below, has a dra- 
matic consequence for the average (exp(i8ips)) . It turns 
out that, while (exp(i5(po^ r )) decays with r exponen- 
tially, the average (exp(j<5<y9j])) falls off only as a power 
law. This, in turn, leads to a slow decay of a zero-bias 
anomaly, 8v(u)/u)\), with u>. 

On the technical level, cancellation of J dxA 2 (x) terms 
leads to a drastic simplification of the disorder averaging 
of Eqs. (|7.10[) , (|7.1ip in the regime II, as compared to the 
averaging of the Friedel oscillations in Section IV Bl since 
the averaging of exp(2i<5(ps) can be performed with the 
help of the Hubbard-Stratonovich transformation. For 
the purpose of functional averaging, it is convenient to 
rewrite Eq. (|9.ip in a slightly different form 



FIG. 14: (Color online) Dimensionless function, T(z), describ- 
ing the shape of a zero-bias anomaly in regime I, is plotted 
from Eq. (|8.6[1 versus dimensionless energy, z = lu/ujo. Inset 
in the lower-right corner: enlarged plot of 1(z) in the domain 
3 < z < 5. 



The steps leading from this expression to the asymptote 
Eq. (|8.9p are outlined in Appendix E. 



IX. ZERO-BIAS ANOMALY IN AVERAGED 
DENSITY OF STATES IN REGIME II 



^(r 2 -n) L Vf 

~\[^- f 2 dxA(x) 
V r 2 Jo 



(9.2) 



Subsequent integration by parts yields the further sim- 
plification of Eq. (|9~2f 



A. Three polarization operators: Averaging of the 
net magnetic phase factor over realizations of 
random magnetic field 



To derive analytical expressions for (uj) and 

5u^~\u) one has to perform averaging of Eqs. (|7.10[) . 
(|7.11[) over realizations of the random field. Such an aver- 
aging has already been carried out for the Friedel oscilla- 
tions. In the latter case we had averaged (exp (2i<5<^ ^r)) ■ 
In the case of the density of states, the exponents to 

be averaged are (exp (2i5(p^^j) , defined by Eqs. (|7.6[) . 
(I7.7|) . Our most important observation is that the net 
phase 5<p^ ' = 5tpo^ ri + Sip ri ^ r2 + 8ip r2 ^Q does not con- 
tain integrals of A 2 (a;), since they cancel out. This can be 
clearly seen from Eq. (|4.17|) . Instead, is expressed 

via integrals of A(x) in the first power as follows 
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dxA(x] 



1 



x ^ dxA(x)) — 
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r% - n 

r 2 x 2 
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(9.1) 
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1 



%k F {r 2 - n) 



* — / dx (r\ — x)h(x, 0) 
V n Jo 



T2 Jo 



dx {r2 — x)h(x, 0) 



(9.3) 



Now the averaging over realizations of h(x, y) can be per- 
formed by a sequence of standard steps outlined below. 



1. Averaging 



Using Eq. (19. ip we rewrite the definition of average 
(exp(2iS<ps)) by introducing the auxiliary integration 
variable, c 
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exp{2iSip^ '}\ = I dc exp i 



s/2 



$ kl /2 y/r 2 - Fl 



— I dx (r\ — x)h(x, 0) 
n Jo 



— / da; (r 2 — x)h(x,0) 

Tl Jo 



(9.4) 



h(x,y) 



where the averaging (. . . )h(x,y) is defined by Eq. (|5. 18|) . Next we use the following integral representation of the 
^-function in Eq. r " 11 



exp{2i(5(^^ '}\ = / dc exp 



f°° dt 

„2\ / al Act , 



,2tt 



exp 



LV r l Jo ®ok F V r 2 - f X 

~ f 2 dx (r 2 - x)h(x,0)~\)\ 



r 2 Jo ®ok F ^r 2 - n 



Now the integration over c can be performed explicitly, yielding 

r0 ° dt 



exp{2i<5<^4 ^} 



,2tt 



e 1 * 2 / 4 / exp<^ -itv^ 



77 T 1 (n -x)/i(x,0) 

n 7o 



*ofc F V r 2 - r l 

ri f r2 dx (r2 — x)h(x, 0) 



»"2 Jo ^o^fa/^ - ri 



h(x,y) 



(9.5) 



(9.6) 



It follows from Eq. (|9.6[) that evaluation of 

;xp{2i<5< ( 5^ '}\ reduces to the Gaussian averaging 

of the exponent of a linear in h(x) functional, which is 
standard 



pr 2 roo 

exp \ -it dx dyh{x,y)f{x)5{y)\) -(9.7) 

J0 J-oo 



h(x,y) 



expj-^- J dxi J dx 2 f{xi)K(xx,Q,x 2 ,Q)f(x 2 )X , 

where K [x\, 0, x 2) 0) is related to the correlator of 
the random field Eq. (|2.ip as follows K(xi,0,x 2 ,0) = 
/iqK(|xi — x 2 \IQ. Subsequent integration over t yields 
the final result 



(9.8) 



\J l + i lo 2 Jo" dx 1 dx 2 f(x 1 )K(x 1 ,0,x 2 ,0)f(x 2 ) 

As seen from Eq. (|9.6[) the function /(x) in Eq. (|9.7p has 
the form 



/-(*) 



V2 



1/2 

(»"2 - x ) 



— (ri — xWri — x) 
ri 



(9.9) 



Averaging of exp|i<5(p^| is performed similarly, and 
also yields Eq. (|9.7[) with /(x) having the form 



/+(*) 



,1/2 . 

$ofc F V 7 ^ 



n + r 2 



ri 



(ri — x)6>(ri — x) 
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n + r 2 



(n + r 2 - x) 



(9.10) 



We emphasize that expression Eq. (|9.8p is general, and is 
valid for arbitrary, /io and £, i.e., in both regimes I and II. 
For the regime I, we had already performed the averaging 
over realizations of the random field. With regard to 
Eq. (|9.8[) . regime I corresponds to replacement of the 
correlator by unity. In regime II, the distances n , r 2 are 
much larger than £. For this reason, in regime II, the 
correlator in Eq. (|9.8[) can be replaced by ^/2~k^S,5{x\ — 
x 2 ), with 7 defined by Eq. (|5.22[) . Then the averages 

(exp{2i,5v3^ ) }) and (exp{2i(5^ +) }) can be expressed in 
terms of dimensionless ratios 



Qi = 



n 



V6r u ' 



Q2 



r 2 



(9.11) 



is defined by 



where the characteristic length, 
Eq. (523). 

Eq. (|9.8p and analogous expression for (exp{2t<fy?£ }) 
are sufficient to perform the averaging over realizations 
of random magnetic field in Eqs. (|7.10p . (|7.1ip . How- 
ever, averaged Eqs. (|7.10p . (|7.1ip contain the real and 
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imaginary parts 

(exp{2* ( 54 ±) }) = uf{g u Q2 ) + i U±{ Bu 02) (9.12) 

of the average exponents, separately. The expressions for 
Uf~ and U2 readily follow after replacing correlator by 
delta-function and performing integrations over x\ and 
x 2 in Eq. 
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(9.14) 



ut = 



y/5y/ei(k + 1 



(9.15) 
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2jg\g\ 



(9.16) 



Final expressions for the contributions (5v~(u>)) and 
(5i/ + (o;)) to the averaged density of states in the second 
regime are obtained by performing integration over f2 in 
Eqs. (j7TTm> and (|7TTTj) and using Eqs. (PT^l) - (j9~TB)) . We 
present this expression in the form similar to Eq. (|8.3[) 



8v ± {uj) 



= V J d 



(9.17) 



where the prefactor T> is defined as 



V 



6 3 / 4 (7rfc F r„) 3 / 2 



(9.18) 



and the dimensionless functions are the following in- 
tegrals over gi , Q2 
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(9.22) 



where z = lo /uj\ is the dimensionless frequency. The new the second regime in a usual way 
energy scale is related to the characteristic length, r„, in 



(9.23) 
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FIG. 15: (Color online) Dimensionless density of states in 
regime II, J{z) — J + (z) + J~{z), is plotted in the units 
of (vqD) from Eqs. J9A9J, (f9720]l . (pHi"]) , (|022)) versus dimen- 
sionless frequency, 2 = uj/uj\, where Mi is defined by Eq. (I9.23P 
and V is defined by Eq. (f9TTH]l . 



The second regime corresponds to long distances, r„ > £, 
travelled by electron. This is reflected in the fact that 
the frequency io\ is smaller than wo _ the characteristic 
frequency for the first regime. Using Eq. (|5.29p . we 
can establish the relation between ojq and u>i, namely, 
uji ~ ujqE 1 / 6 , where e is the small parameter, defined by 
Eq. (|3.10|) . We emphasize that the second regime exists 
only if the condition e <C 1 is met. 

It is important to compare the scale wi to the "diffu- 
sive" energy scale uj Ai{{ ~ v F /l tI , where l tt is the transport 
mean free path. In the regime II we have 2 ^ 



Mil 



(9.24) 



In this estimate the combination, /iq£ 3 v f /$q, stands for a 
single-particle scattering rate, calculated from the golden 
rule, with /iq^V^o coming from the square of the matrix 
element; the factor (fc F £) accounts for the small-angle 
scattering. Eq. (|9.24p leads to the following relation be- 
tween the transport mes.il free path and. Vn 



(9.25) 



As follows from Eq. (|9.25[) . the distance r„, over which 
the phase of the Friedel oscillations is preserved, is in- 
termediate between l tI and £. Indeed, the ratio l tI /r u 
is ~ k F r u ~ fc P £/-^/i. This ratio is large both because 
k F £ ^> 1 and because e«l. Thus we conclude that the 
energy scale, ui, is much larger than 

^diffi since c^diff 
is ~ r n /l tT <C 1, i.e., the conventional diffusive zero-bias 
anomaly develops at frequencies much smaller than the 
width of the zero-bias anomaly in regime II. 



B. Discussion 

Dimensionless density of states, J = 3~ + J + , is plot- 
ted in Fig. [T51 It is seen that the function J(z) exhibits 
pronounced minimum at z ss 0.75, which is followed by a 
monotonous decay. This behavior should be contrasted 
to the dimensionless density of states in the regime I, 
plotted in Fig. [14] The difference is that the function 
I exhibits damped oscillations with alternating maxima 
and minima, while J contains only a single minimum. 
This difference is not unexpected on qualitative grounds. 



Indeed, the distance 



at which the oscillations are 



formed in regime I, is much smaller than the correla- 
tion radius, £, while the characteristic distance, ~ r„, in 
regime II is much bigger than £. Therefore, it is remark- 
able that J7(z) exhibits even a single minimum. However, 
qualitative difference between T(z) and J{z) at large z 
is much harder to trace from their respective representa- 
tions as double integrals over p\ and p2 [see Eqs. (18.61) . 
ipn$|) , (j9~2T]|) . dOTJ, (gUl) ]. The structure of one of 
several contributions to T(z) and J(z) can be loosely 
rewritten as 
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The integrands in Eq. (|9~2^|) and Eq. (f9~2"7) differ only 
by the structure of the denominators. This difference 
can be traced to Eq. (|9.8p in which the correlator is set 
either constant (regime I) or a 5- function (regime II). 
From the form of the contribution Eq. (|9.26[) . it is not 
obvious at all that the large-z behavior is determined by 
well-defined values pi = p 2 = po in the complex plane, 
with po satisfying 1 + Pq = 0, so that the contribution is 
oscillatory Eq. (|8.9[) . This fact was established above by 
taking the large-z asymptote prior to the averaging over 
realizations. It is also supported by numerics in Fig. [TJ] 
Monotonous behavior of 3{z) at large z implies that 
the integral Eq. (|9.27[) is not dominated by distinct com- 
plex pi = p 2 = po, such that 1 + p~Q = 0. The 
only vague explanation of this is that denominator, 
+ plpKpi + p 2 ) 2 , in Eq. (19.26[) fixes pi ~ p 2 ~ Po 
much more efficiently that the denominator, yl + p\p\, 
in Eq. (|9.27[) fixes pi, p 2 near po- 



X. IMPLICATIONS 

A. Half-filled Landau level 

Experimental situation of a two-dimensional electron 
gas placed in inhomogencous magnetic field can be cre- 
ated artificially, see, e.g., Refs. 11(0130000111 



23 



Il8l . This situation also emerges in electron gas in a 
strong constant magnetic field, when the filling factor of 
the lowest Landau level is close to 1/2. In the latter 
case, constant field transforms electrons into composite 
fermions 8 -^, with well defined Fermi surfac o 67 i 68 ' 69 i 70 ' 71 , 
while the randomness of magnetic field is due to spa- 
tial inhomogeneity of the electron density. Transport 
properties of noninteracting gas of composite fermions 
under these conditions were considered theoretically in 
Refs. [lllllllllllllllllS. 

With regard to the tunnel density of states near the 
half-filling, for the case of ho mogen eous gas, it was ad- 
dressed theoretically in Refs. I72ll73ll74l both for tunnel- 
ing into the bulk and into the edge. Unlike interact- 
ing homogeneous electron gas* 4 - composite fermions are 
expected to exhibit a zero-bias anomaly even without 
inhomogeneit y 72 ^ 73 ^ 7,1 . This difference between composite 
fermions and free electrons can be traced to the form of 
density-density correlator of composite fermions at small 
momenta^. Namely, the pole of this correlator defines 
the mode of neutral excitations with dispersion u> tx iq 3 , 
even slower than the diffusive mode in the presence of 
disorder. Resulting suppression of tunneling into the 
edge of homogeneous electron gas at half filling, predicted 
in Refs. [73], [3, turned out to be stronger than in the 
experimen t 75 ' 76 . 

It is convenient to express random static magnetic field 
originating from spatial inhomogeneity with magnitude 
Sn, in the units of the cyclotron frequency 



Suj c 

n 1 / 2 



2Sn 



(10.1) 



where n\/2 is concentration of electrons at which the 
filling factor in the field, Qi/ 2 , is equal to 1/2. Den- 
sity fluctuations not only smear out the "intrinsic" zero- 
bias anomaly, but also give rise to the smooth-disorder- 
induced zero-bias anomaly, studied in the present paper. 
Quantitatively, we predict the following relation between 
the width of zero-bias anomaly and the magnitude, 5n of 
the density fluctuations 
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This relation follows directly from Eq. (|3.11j) and applies 
for smooth fluctuations with spatial scale, £, satisfying 
the condition 
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This condition is equivalent to the condition e > 1, where 
the parameter e is defined by Eq. (|3.10J) . In the opposite 
case of "fast" fluctuations the width, u\, is given by 



^1/2 [«-l/2£ 2 ] 
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as follows from Eq. (|3 . llj) . Concerning the magni- 
tude of the anomaly, Eqs. (|8.5[) and (|9.18[) predict 
b~vjvQ ~ (Sn/ni/ 2 ) for slow fluctuations Eq. (|10.2p . and 



.3/2 



ni/2? 2 



3/8 



for the fast fluctua- 



6v/i>o — (Sn/ny 2 ) 
tions Eq. (|10.4|) . respectively. 

Qualitative difference between the "intrinsic" zero- 
bias anomal y 72 ' 73 ' 74 and inhomogeneity-induced zero- 
bias anomaly, considered in the present paper, is that 
the latter necessarily involves electron-electron scatter- 
ing processes with momentum transfer w 2fc F . As was 
mentioned above, the intrinsic anomaly gets stronger to- 
wards the edg o 73 ' 74 . We would like to emphasize that 
the anomaly due to the 2&; F -processes also gets stronger 
towards the edge. The reason is that the average electron 
concentration decreases monotonically upon approaching 
the edge. This decrease translates into a non-fluctuating 
magnetic field, acting on composite fermions^ 3 -, which in- 
creases towards the edge. Correction, 8v(u>), to the den- 
sity of states in this case is given by Eq. (|7.13p . and is 
plotted in Fig. [101 Then we conclude that the ratio of 
magnitudes, 5u hvlv /Su^, is simply ~ {5n bulk /Sn odso ) < 
1, where <5n bulk and Sn odgo are the deviations of electron 
density from ni/2 in the bulk and at the edge, respec- 
tively. The widths of 5v hum (u>) and <5^ odgo (a->) are related 

as — (<focd g c/c>n-buik) 2/3 < I- 



B. Spin-fermion model 

Similarly to composite fermions, the dispersion of 
neutral excitations right at the critical point in the 
spin-fermion model is dominated by a slow mode , 77 ' 78 
u> oc iq 3 . Outside the critical region, the propagator of 
the neutral excitations (bosons) in the spin-disordered 
phase has a conventional Ornstein-Zernike form xil) ^ 
1/ (q 2 +£,~ 2 ), where £ is the correlation radius, which 
diverges at the critical point. Interaction of electrons 
with slow critical fluctuations can be viewed as scatter- 
ing by the smooth disorder. The question that we will 
discuss below is how the growth of £, upon approaching 
the critical point, manifests itself in the behavior of the 
averaged (over the fluctuations of the order parameter) 
density of states. Our calculations demonstrate that the 
dimensionless parameter e, defined by Eq. ()3.10j) . plays 
a crucial role. 

Traditionally, in the studies of the response func- 
tions, like spin susceptibility, of two-dimensional elec- 
trons near the quantum critical point, see, e.g., Refs. 
77J8 ; 79.80 ; 81.82, electrons are treated as ballistic. More 
specifically, they interact only with critical fluctuations, 
but not with each other. Transport at the quantum 
critical point was also considered for non-interacting 
baUistioS 3 - or diffusive^ 4 - electrons that are scattered by 
bosonic excitations. 

In all theoretical treatments of the spin-fermion model, 
modification of the response of the electron gas due to in- 
teraction with bosons was governed by the processes with 
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small momentum transfer. Our main point is that in- 
corporating direct electron-electron interactions into the 
spin-fermion model gives rise to a novel feature in the re- 
sponse of the electron gas near the critical point in spin- 
disordered phase. The underlying reason is that, while 
critical bosonic fluctuations are "smooth", so that their 
momenta are <C k F , electron-electron interactions allow 
2fc F -processes. Then the physics, discussed in the present 
paper, emerges in the following way: 

(i) interaction with slow bosonic fluctuations, curves 
slightly the electron trajectories; 

(ii) interaction between the electrons, moving along 
slightly curved trajectories, generates a small energy 
scale, which reflects the "degree" of curving; 

(iii) the degree of curving grows with correlation ra- 
dius, £, of the bosonic excitations. 

As a result, the character of critical fluctuations is re- 
flected in the density of states, 5v(pS), in a very non- 
trivial fashion. Namely, they give rise to the lively 
low-frequency feature and even aperiodic oscillations in 
5v{ui), as was demonstrated above. This suggests that 
information about proximity to the critical point can be 
inferred from tunneling experiments. 

To quantify the above scenario, we will assume for 
simplicity^ that bosonic critical fluctuations of magne- 
tization, S(r), interact with electron spins not as <x ■ S, 
where cr are the Pauli matrices, but via the position- 
dependent Zeeman energy, E z (r), with characteristic 
magnitude, Eq. Assuming that the fluctuations, S(r), 
are static, we get for correlator of random Zeeman en- 
ergy, E z (r), the standard expression 



(E z (r)E z (r>)) = E 2 J g 



q 2 + t~ 2 
= El K - x 2 |/0 , (10.5) 

where Ko is the Macdonald function. 

As a next step, we notice that the force, \7E z (r), curves 
the electron trajectories in the same way as random mag- 
netic field, h(x,y). This allows us to use general expres- 
sions Eqs. (|7.10[) . (|7.1ip for the interaction correction 
to the density of states. We can also employ the re- 
sult Eq. (|9.8[) for the general averaging procedure, i.e., to 
treat critical fluctuations as a disorder. With the help of 
Eq. (|10.5|) the result Eq. (|9.8[) assumes the form 



1 + iE 2 



exp{2i(5</52 ^} 
x f-{xi)f-(x 2 )d Xl d X2 K 



'o jo 
\xi - x 2 \ 



dx\dx 2 (10.6) 



-1/2 



where the function /_ is defined by Eq. (|9.9p for the 
case of random magnetic field. For the case of random 

1 /2 

Zeeman energy, the prefactor, I/^o^f , should be re- 

1/2 

placed by k F /E F . Characteristic energy scales can be 
now inferred from Eq. (|10.6[) on the basis of the following 
reasoning. Characteristic distances r\, r 2 in Eq. (|10.6|) 



are determined by the condition 



i, (10.7) 

^0 



where we performed integration by parts in Eq. (|10.6[) . 
Then the characteristic width of a zero-bias anomaly is 
equal tow~ Up/ri ~ v F /r 2 . 

Recall now, that in the case of random magnetic field, 
double integral in the left-hand side of Eq. (|10.7p did 
not contain derivatives and was oc r\ in the regime I, 
and oc r|£ in regime II, respectively. This is because the 



function, f_(x{), is ■ 



„l/2 



at X\ ~ r 2 , see Eq. (|9.9[) . Due 



to the fact that the effective "force" in the spin-fermion 
model is oc V£' z (r), the left-hand side in Eq. (|10.7|) is 
~ k F r 2 /E 2 for £ > r 2 . In this limit, Eq. pTJTTjl yields 
(with logarithmic in ^fr 2 accuracy) 



E v 



(10.8) 



Note that E c is independent of £. We conclude that, upon 
approaching the critical point, as the correlation radius 
exceeds the value £ c , the zero-bias anomaly "freezes". 
Its form is shown in Fig. 1141 and its magnitude is 
~ (Eq/E f ) 3 . An alternative way to recover the scales 
Eq. (|10.8p is to notice that parameter e, which is defined 
by Eq. (|3.10p in context of random magnetic field, in the 
situation with random Zeeman energy acquires the form 
e = (k F £){E /E F f. Then £ c given by Eq. (flTl8)) cor- 
responds to e = 1, i.e., to the boundary of the regime 
I. 

For £ < £ c the integral in the left-hand side of 
Eq. (|10.7p is proportional to £ and is independent of r 2 . 
Then Eq. (|10.7p does not have a solution. Therefore, 
characteristic r\ and r 2 in the expression for the den- 
sity of states are ~ £, and the width of the anomaly 
is simply ~ v F /£ — E c (t; c /t;). Concerning the magni- 
tude of the anomaly at £ < £ c , it should be estimated 
with the account that the integral in right-hand side 
of Eq. (|10.6p is smaller than 1 for all r 2 . Therefore, 

(exp{2i5(p^ ^}) in Eq. (|10.6f> can be approximately re- 

placed by {l-(i/2^) [r 2 9(£ - r a ) + £6(r 2 -£)]}, where 
the second term is a small correction. However, only this 
correction causes a zero-bias anomaly. Substituting this 
correction into Eq. (|7.11[) . we find the estimate for the 
magnitude, 



5v 



Eo 
E F 



1/2 



E F 



(fc F £)V2 



(10.9) 



We conclude that, as £ grows and approaches £ c , the mag- 
nitude of the anomaly falls off as 1/ v^, and the anomaly 
narrows as l/£. 

The remaining issue to discuss is whether or not the 
assumption that fluctuating Zeeman energy, E z (r), is 



25 



static applies at relevant frequency and spatial scales, 
E c and £ c . For this purpose, we recall the correla- 
tor of Zeeman energies in the momentum space does 
not have a simple Ornstein-Zcrnikc form but is rather 
(|i? z (q)| 2 ) cx 1/ (q 2 + £~ 2 + <;uj/q), where the dynamic 
term, suj/q, describes the damping of bosons due to 
creation of electron-hole pairs. The prefactor q (the 
Landau damping coefficient) is thus quadratic in cou- 
pling of electrons to the spin density fluctuations, i.e., 
<r cx Eq. For characteristic frequencies the dynamic 
term, q(u)/q) ~ <;E C £ C ~ c;v F . Therefore, it is negligible 
only if the condition, = k 2 , (En/E F ) 4 3> SV F , holds. 
With <; being proportional to Eq, the above condition 
is met for large enough coupling, Eq. In the opposite 
case, when the dynamic part of correlator dominates at 
u) r~j E c and q ~ , the zero-bias anomaly develops 
only away from the critical point when £ becomes smaller 
than (<^w F ) -1 / 2 . Upon further departure from the critical 
point, our prediction Sv/vq cx £~ 1 / 2 and ui ~ v F /£ should 
apply. Note finally, that, directly at the critical point, the 
slow mode to ~ iq 3 gives rise to the "intrinsic" zero-bias 
anomaly^ similar to the composite fermions. 
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APPENDIX A: POLARIZATION OPERATOR IN 
THE COORDINATE SPACE 



Eq. (|A2|) with respect to frequency 

0(qv F — ui} iQ(u> — qv F ) 



\Jq 2 ~ {uj/v f ) 2 y/ (uj/v f ) 2 - q 2 



ds Jo(qs) exp 



(A3) 



Substituting Eq. (|A2[) into Eq. (|A3|) and using the or- 
thogonality relation L dq qJo(qs)Jo(qr) — S(r — s)/r, 
we readily obtain 



n r,u;) = — exp }. (A4) 

v F r v F 



In order to calculate Il2fc F (r, ui) we use the form of polar- 
ization operator in momentum space for \q — 2k F \ <C k F 
and ui <C E F 



n2/t F (q, w) = 



i 



(A5) 



\sj q — 2k F + lu/v f + \j q — 2k F - u)/v F ^j 



where the square roots should be understood as \fx — > 
sign(x)v / a;. Then the integral over q in Eq. (|A1|) assumes 
the form 



n 2feF (r,w) = -v I dqqJ (qr) ^ q - 2k F + uj/v f 
Jo L 



+yj q — 2k F — uj/v f 



/4fe F / ~ / 7i- 

/ dq cos or 

tit In \ 4 



\J q — 2k F + uj/v f + \fq~— 2k F — tu/i 



(A6) 



Here we derive Eqs. (|4.2I) and (|4.3[) for polarization op- 
erator in coordinate space using the known expression 59 
for n(q, uj) in the momentum space. Since we are inter- 
ested in behavior of U(r,u>) at distances |r| ^> fc^ 1 , it is 
sufficient to perform the Fourier transform 



U(t,u) = ± Jdqe^ n(q,«), 



(Al) 



using the asymptotes of II(q, ui) at small q <C k F and at 
q close to 2fc F . The small-g asymptote of II(q, ui) has the 
form 



where we used that fact that k F r ^> 1 and replaced the 
Bessel function by its large-g asymptotics. Integration 
over variable q in Eq. (|A6|) is performed with the use of 
the identity 



dz cos z\[z — a = —— sin (^a + — ^ , (A7) 



and yields the zero-temperature limit of Eq. (|4.3 



APPENDIX B: POLARIZATION OPERATOR IN 
A CONSTANT MAGNETIC FIELD 



n (q,w) 



iu Q(qv F — ui) 
y/ q 2 v 2 - ui 2 

to Q(uj — qv F ) 
y/ui 2 - q 2 v 2 



(A2) 



where Q(x) is the step- function. The easiest way to per- 
form the integration Eq. (|A1[) is to first Fourier transform 



We start from the general expression^ for the polariz- 
ability in arbitrary magnetic field 



n(?) 



~ OO OO 

Zra sr^ sr^ 



ni—O U2—0 



exp(-a 2 Z 2 /2) (^fj (ff) , (Bl) 
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where L™^ ni (x) and L™* ™ 2 (a;) are the Laguerre poly- 
nomials, and /„ = |exp[(n — N F )fko c /T\ + 1} 1 is the 
Fermi distribution. At small q <C k F Eq. (|B1[) yields 63 
II (q) = — (m/n)[l — J^qR^)], i.e., the characteristic 
scale is q ~ -R^ 1 - For (g — 2fc F ) <C & F it is convenient 
to perform the summation over the Landau levels with 
the help of the following integral representation of the 
Laguerre polynomial 



2tt 



271 



(W 



o (1 - e M ) 



re 



it) 



— imO 



(B2) 



In the vicinity q = 2k F Eq. (|B2[) contains a small fac- 
tor exp(— q 2 l 2 /2). This factor is compensated by the 
product of Laguerre polynomials, since each of them is 
oc exp(x/2), which comes from the exponent in Eq. (|B2|) 
taken at 9 — ir. With contribution from the vicinity 9 = 
7r dominating the integral (|B2p . we can expand the inte- 
grand around this point as exp[x/2+iirm+i(f>(ip)] /2 n+1 , 
where if) = [B — ir), and the phase, 4>(ip), is equal to 



4>W 



r n + 

m 

4 2 



rip 3 
~48 



(B3) 



Now we make use of the fact that only relatively small 
number ~ (fc F Z) 2 / 3 <C N F of Landau levels around E F 
contribute to the sum Eq. (|B2p . This suggests that we 
can present n\ and n 2 as n\ = N F +mi and n 2 = N F — m 2 , 
respectively, and extend the sum over mi, from —00 
to +00. After that the summation over Landau levels 
can be easily carried out with the help of the following 
identity 



E 



fN F — mi fN F +ni2 



mi ,m2 — — 00 



x cos [(mi — 7712)0: + 0\ — 



mi + m 2 

2tt 2 Tcos/3 
Hue sinh(27r|a|T/fia' c ) 



(B4) 



As a next step, we substitute the representation Eq. (|B2|) 
of Laguerre polynomials with integrand expanded ac- 
cording to Eq. (fB3j) . into Eq. (|Bl|) . Upon this substi- 
tution, we perform the summation over Landau levels 
using the relation Eq. (|B4jl . Then the double integral, 
which emerges in Eq. (|B 1 1) as a result of representing the 
two laguerre polynomials Eq. (|B2[) . assumes the form 



dipidip 2 



x cos 



-00 1^1 +'02 I 

' I) ±- 

12 



(V>i + V>: 



+ ^2) 



(B5) 



where Sq = q — 2k F . Note, that integration over the differ- 
ence, (tpi — 1P2), in Eq. (|B5|) can be performed explicitly. 
It is convenient to present the final result not for H(q), 
but rather for the derivative, W(q,T) = dU(q,T)/dq. 
Knowledge of II'(g,T) is sufficient for finding the large- 
distance behavior of the potential, created by the short- 
range impurity. Indeed, this potential can be expressed 
directly through II'(2fe p + Q) as follows 



V H (r) 



V(2k F )g 
2(nk F r) 3 / 2 



dQ sin 



(2fc F - 

n'(2fc F 



,T). (B6) 



At zero temperature and in a zero magnetic field we have 
II'(g,0) cx B{8q)/\fSq. At finite magnetic field and finate 
temperature, taking derivative of Eq. (|B5|) with respect 
to Sq, we arrive to the result 



n'(«,r) 



dr x 1 / 2 



2 1 / 3 mT 

(7rfc F p ) 1/2 £o Jo sinh(27rxT/e ) 



s 5q 



1 



. SUI 1 2 2 / 3 — x + -x u + , 
Po 3 4 



(B7) 



In the limit T -> 0, substitution of Eq. (B7j) into Eq. JBSJ 
and integration over Q reproduces Eq. (|3.ip . 

Interestingly, for T = 0, the integral Eq. (|B7|) can be 
evaluated analytically 



n'(«) 



Ai 



Po/ VPo/ 



(B8) 



where is the Airy function, and Bi(z) is another 

solution of the Airy equation defined, e.g., in Ref. l86l . 
It is seen that the singularity at q = 2fc F is smeared by 
the magnetic field in a rather peculiar way: for positive 
Sq ^ Po the (Sq) -1 / 2 zero- field behavior [see Eq. (|A6j) ] is 
restored. However, for large negative Sq/po, the deriva- 
tive II' (g) approaches zero with oscillations, namely, as 
cos[4(| < 5g|/po) 3/2 /3]/(|(5g|) 1 / 2 . As the difference 2k F - q 
increases and becomes comparable to k F , these oscilla- 
tions cross over to the "classical" oscillations^ 3 - II' (q) oc 
J (qR L )J 1 (qR L ) cx cos(2qR L ). 



APPENDIX C: EVALUATION OF THE 
FUNCTIONAL INTEGRAL 



Upon combining Eqs. (|5.24[) and (|5.26p the quadratic 
form in the exponent in the numerator of the functional 
integral Eq. (|5. 18|) assumes the form 
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2iScp(r) - W{h} = 



2ier 6 1 



C 3 12 



I dqA^.q + y^c ra | I dqAn, q + f dq A a . q G{A n } 

J n>0 ^ ^ 



(CI) 



2r ^ /■ |A, 12 



n>0 



^ ^ /,/„ 2 



JC(g) 7U /C(«) 



with numerical coefficients c„ = l/2-7r 2 n 2 and & n = — c„ + «/27rn denned by Eq. Q5.27[) . In the above expression we 
had introduced a short-hand notation 



Q{A„, q } = J2[ bn dq An 'i + b *n dq -K 

n>0 J J 



(C2) 



We adopt the following sequence of integration over the variables A n , q - First we integrate over Ao, q using the following 
decoupling 



H{<?} = J [J <Mo, t cxpj I dq + |J y dg A« 2 +iGjdq A , 



(C3) 



3r 3 



2 ^3 7 n^,^oexp|-5||o+iBo/^A,-^/^ 



' iG / d,y.4,,„ 



r 



where we had introduced an auxiliary variable Bq. Func- where c„ is related to c„ via a dimensionless factor 
tion H{£/} combines all integrals in Eq. (|C1[) containing 
Ao iq . Subsequent integration first over the variables Ao tQ 
and then over the auxiliary variable Bq yields 



H{6} = 



/ 7rf 7 exp{iF{r)g 2 } 
irJdqIC(q) _ ^ ^7^' 



where we had used the definition r n = 2£/ (-\/27r7e) 
In Eq. (|C4|) the complex function .F(r) is defined as 



^(r) 



3 £ V2 Ujdqt{q) 



3/2 
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(*)' 



24i — 



(C4) 



As a result of integration over Ao, q the exponent in the 
functional integral Eq. (|C1[) assumes the form 



^ cj / A n ,, 

~. ^ n 



n>0 

2r 



n>0 



2er 3 



(C6) 



The first and the third terms in Eq. (|C5[) contain squares 
(C5) °f ^ ne nnear combinations of A n . q - To decouple these 
squares, we introduce a set of auxiliary variables, a n ,a* 
for the first term, and one auxiliary variable, «o for the 
third term as follows 



J 



iV\fdqA n , q \ 2 = 



i Jda n da n exp{-z|a„| 2 + V 1/2 a* ^ .4„,« - V 1/2 a n J dq A* hq j, 



(C7) 
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^(r)S 2 _ 



dao exp 



ij4 
AT{r) 



(C8) 



Note that Im [l/.F(r)] < 0, so that the decoupling Eq. (|C8j) of the quadratic in Q term in the exponent of Eq. (|C4|) is 
justified. 

As a next step, we perform gaussian integration over the infinite set of variables, {-4„ !g } 



dA n , q dA* nq exp{ I dq 



2r\A n , q \ 2 
' 7^(<?) 



(cl/ 2 a* n + ia b n ^ +A* nq (^~clJ 2 a n + ia b 



2l7T 



-exp|- 



icj a* + a b r 



(2r/ 7 0fdq [lC(q) 

As follows from Eqs. (|C7[) and (|C9[) . the integrals over all a n are gaussian and can be easily evaluated 

2 7? 



da n da* n 
2V 



■ exp 



j-i|a„| z 



ic l J 2 a* n + a b r 



2r 



2r 



2r - ij^Cn J dqfC(q) 



expi 



7-foo 1 K 
2r 



dqit(q)- 



dq IC(q)\ 

1 2 

Ur 2 +2n£c n JdqK,(q)\ 



(C9) 



(CIO) 



The remaining integral over ao is also gaussian. Note now, that the denominator in Eq. (|5.18p . responsible for the 
normalization, can be evaluated by performing the same steps as above. This evaluation amounts to setting c n = 
in Eq. (|C10|) and taking the limit r„ — * oo in Eq. (|C4[) . As a result, the functional integral reduces to the ratio of the 
ordinary integrals 



(*)* 



If n - 

AA n 2 -2i(r/r n ) 2 /n 



f^jlao exp {-uia - Mia 2 ,} 
S-oo da o ex P {-wa - u al} 



(Cll) 



where the coefficients w, and uq are defined as 



"o 



2r 
r 



n>0 

dqlC{q) 



dqfC{q), 

3/2 



while the definition of the coefficient u\ is the following 



Ui 



AT{r) 



E/ M 



7^ / dqK{q) 
c„\b n \ 2 



/ 4i r 2 -)- 2r7^c„ / dqJC(q) 



(C13) 



For characteristic r ~ r„ the first term in Eq. (|C13[) is 



-1/2 



as follows from Eq. (|C4[) . On the other hand, 



Next we notice that, for r ~ r„ both uq and ui are of 



the same order and are 



-1/2 



On the other hand, 



as seen from Eq. (|C12[) . the parameter w for r ~ r u is 
small, ui ~ e 1 / 2 . This allows to disregard w both in 
numerator and denominator in Eq. (|C11[) . so that the 

(C12) ratio of integrals reduces to (uq/ui) 1 ^ 2 . Using Eq. (|C4[) . 

this ratio can be rewritten as [l — (2i/3)(r/r„) 2 ] . 
Substituting it into Eq. (|C11|) . we arrive at Eq. (|5.28|) in 
Section N\ 



the product r£c n in the denominator of the second term 
in Eq. (|UT3|) is - er 4 /£ 2 ~ r 4 /r 2 . Thus, for r ~ r„ 
both terms in the denominator of the sum in the second 
term are ~ r 2 . The numerator in the sum over n is 
~ £~ x / 2 for r ~ r„. Then the estimate for the second 
term in Eq. jCT3|) is £, 2 /r 2 s 1/2 , so that the second term 

is smaller than the first term in parameter £ /r 2 ~ e. can be naturally divided into two parts 1(z) = T + +I_ 



APPENDIX D: ANALYSIS OF THE INTEGRALS 
EQ. ([OD 



The dimensionless function T{z) defined by Eq. 
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where 



T_(z)=/ f^^a/s I dz' ain[(z ~ z')(p 2 + Pl )] 

Jp 2 > Pl (P2Pl) 3/2 Jo 

x VP2 - Pi I sin[7r/4+ (z + z')(p 2 - pi)] 



l + ^l + P lp\{p 2 -pi) 



i + pIp\{p2 - pi) 2 



+ cos[tt/4+ (z + z')(p 2 -/n)] 



V 1 + pIpKpz - pi) 2 - 1 



l+P|Pl(P2 ~ Pi) 



2 ( ' 



(Dl) 



and 



l + (z)=[ T^-rin [ dz'sm[(z- z')(p 2 + Pl )] 

Jp2> P l (P2Pl) < J0 

Xy/p2 + Pi\ sin[7r/4 -(z + z')(p 2 + pi)] (D2) 



1+^1+ P jp((p2+Pl) 2 



1 + P 2 2 P 2 i(P2 + Pi) 2 



+ cos[tt/4- (z + z')( P 2 +pi)] 



/ yi + ^(P2 + pi) 2 -i l 
V ^+pIp\{p2+pi? y 

The complexity in numerical evaluation of X + and I_ 
stems from the fact that, upon integration over z' , both 
integrals turn into the sums of two contributions, each 
of which is divergent in the limit z — > 0. Therefore, it 
is necessary to rewrite the result of integration over z' 
in X + and in X_ in such a way that cancellation of the 
divergent contributions is explicit. 

We start with Z_. Integration over z 1 generates the 
combination of three terms 



Pi + P2 

P1P2 
1 

P2 



■ COS 



- + 2z( Pl - p 2 ) 



— cos + 2zpx] — cos - 22^2) (D3) 

p 2 V4 'pi ^4 / 



In order to treat all these three terms on the equal foot- 
ing, in the first term of Eq. (|D3|I we introduce the follow- 
ing new variables 



z = z(p 2 - pi), 
P2P1 



X 



Z 3 (p 2 ~ Pi) 2 



(D4) 



In the second term we introduce z — zpi, and finally, in 
the third term, z = zp 2 . After that, the expression for 



X_ assumes the form 
1 



l.(z) 



dx 



V2z 3 J x 5 / 2 



1 



dx 



V2z 3 Jo x 5/2 

1 r 1 ^ 3 dx 



fi(as)-Fi(0) 
F 2 (x,z) (D5) 

F 3 (x,z)-F 3 (0,Q) , 



Vlz^Jo ^ 5/2 
where the functions F\, F 2 , and F3 are defined as 



Fi(x) = j -^rj^{ (cos22-sin2SV 



^5/2 I *~~) U 1 _|_ p x 2 



(cos 2z 



+ sin 2z) \ - 



1 + Vi + z 6 a; 2 



+ z 6 z 2 



(D6) 



^1 + ^3-1 




f xz 3 


r^LJ 

7o ^ 5/2 l 


(sin 25 



(D7) 



cos zz) 
(cos2z + sin2z) 



Wl + S e x 2 + 1 
1 + z 6 a; 2 



'vTTIV-l 



1 + z b x 



56^,2 ( ' 



F 3 (a;,z) = - 



1 3 1 



2j|-^ 3 



dz 

1^2 ' 



/ „. . „_s /VT + 5 6 a; 2 + 1 
(cos 2z + sm 2z) W 1 - |^ 



+ (cos2z-sin2^ V \ + + ^ ^ (DM 



Subtraction of x ~ values from F±(x) and i<3(2;,z) in 
Eq. (|D5|) insures the convergence of integrals over 2 in 
Eqs. (|D6|) and (|D8|) . On the other hand, this subtraction 
shifts Z_ by z-independent constant. 

It is seen that, in the limit z — > 0, the difference 
F 2 (x,z) — F2(0,z) behaves as z 3 , so that the contribu- 
tion from F 2 to T-(z) remains finite in this limit. On 
the other hand, the contributions from F\ and F3 both 
behave as 1/z 3 . To demonstrate that the two divergent 
contributions cancel out, we divide the integration do- 
main in the first term of X_ into the intervals {0, l/4z 3 } 
and {l/4z 3 , 00}. We then combine the two integrals from 
to 1 /4z 3 to obtain 
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and 



2L(*) 



^ dx 



+ (f 3 (x,z)-F 3 (0,z)) 
x Fi(a;)-Fi(0) 



,5/2 



-.Fi(0) 



1 



°° dr 

F 2 (x,z). 



V2z 3 Jo x 5 / 2 



(D9) 



r(2) 



2 r 1 '^ 3 
~ V2J0 

(sin 25 — cos 2z) 



(cos 25 + sin 25) 



a;3/2 A /i _ xz s J 



dz, 
J3J2 



Vl + z 6 x 2 - 1 

1 + Z 6 X 2 



vm w x 2 + 1 



-V2 



(D13) 



1 + 5 6 .t 2 

Both these contributions arc finite in the limit z — > 0. 



The second and the third terms in Eq. (JD9J) are conver- 
gent in the limit z — > 0. The integrand in the first term 
has the form 



Fx{x) - F(0) + F 3 (x, z) - F 3 (0, z) 

3 



(D10) 



= <l- 



fiz 



(cos 25 + sin2i 



1 + Z 6 X 2 



V2 



(cos 25 — sin 25) 



' VI + z 6 x 2 - 1 
1 + z 6 x 2 



We see that in the limit z — ► expression in the curly 
brackets behaves as oc z 3 , and thus cancels the divergent 
prefactor. Now all three terms in Eq. (|D9[) yield a finite 
contribution at z — > 0. Our numerical results for T-{z) 
were obtained from Eq. (|D9|I . 

We now turn to Z+(z). In order to deal with small- z 
behavior in the integral Eq. (|D2|) . we introduce, after per- 
forming integration over z' , the following new variables 



5 = z(p 2 + pi), 
P2P1 



X 



Z 3 ( P 2+Pl) 2 



(Dll) 



Then one obtains T+ = T\ + T\ , where the two contri- 
butions are given by 



r(l) 



1 

V2 Jo 



1/42 



dx 



x 3 l 2 \l \ — 



dz 



3 Jo Z 



5/2 



(D12) 



(l — sin 25 — cos 2Z 



. 1 — cos 25 + sin 25) 



Vl + 5 6 x 2 - 1 
1 + S 6 x 2 



Wl + S 6 x 2 + 1 
1 + z 6 x 2 



V2 



APPENDIX E: ANAL YSIS O F TH E INTEGRALS 
EQS. ([7723)1 . (f7T24|) 



In the integral Eq. (|7.23|) we perform the following 
change of variables 



v + A 



after which it acquires the form 

3 . 2 13/6 ,ao dy 



tt3/2 / ,,1/2 / r 3/ 



(El) 



cos 



7T Z 

iz 

4 v + 4 





7T" 


— COS 


a;z — - 




L 4J 



In the integral Eq. ()7.24|) we perform the following 
change of variables 



"=2 1 



P2=2 



u + 4 



v + 4 



after which it acquires the form 
P 1 -(x) = 



r 3/2 



dv 


,00 ( 


( v + 4)1/2 y o z 




7T" 


— cos 


iz + - 




4J 



dz 



(E3) 



(E4) 



It is convenient to present / °° dv in Eq. (|E4|) as the fol- 
lowing difference of integrals 

3 ■ 2 13 /6 i-oo dz 

P l ( X ) = =372 / — ( E5 ) 



dv 



dv 



Vv + 4 J_ 4 V« + 4 



x cos 



z 3 ~ 






7T" 




— COS 


xz - 




1> 






" 4. 
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We now observe that the the second term cancels iden- 
tically the function P + . Then we readily arrive to 
Eq. fT2H]l . 



APPENDIX F: ASYMPTOTICS OF THE 
DENSITY OF STATES 



The argument of cosine in Eq. (|F2[) has a sharp mini- 
mum at p = po = (2 13/6 /3 1/2 ) y/u/uh, which allows to 
perform the integration over p by introducing dp = p — po 
and extending the integration over Sp from minus to plus 
infinity. This yields the following asymptote of 8v{u>) 



The idea of derivation of Eq. (|8~9|) from Eq. (pHO)) is 
that the major contribution to the integral Eq. (|8. 10|) 
comes from the domain |p 2 — Px\ Pi,P2, i-e., from the 
domain where p\ and P2 are close to each other. To 
make use of this simplification we rewrite the argument 
of cosine in Eq. (|8.10[) as 



(Pi +P2? 



TV 

4 



2 7/3 



UJ 



{pi + P2){P2 - Plf 



(Fl) 



where we had introduced tOh = wo(/i//io) 2 ^ 3 - It is seen 
from Eq. (|F1|) that the typical value of (p2 + Pi) is 
(u / ujh) 1 ^ 2 3> 1, while the typical value of (p2 — Pi) is 



(P2 + Pi)~ 1/2 ~ (w/wh) -1 ' 4 , i.e., the relevant difference 
P2 — Pi is small indeed. This allows to extend the inte- 
gration over p2 — pi from zero to infinity and perform the 
integral. This yields 



5v{lo) 



(voVfuto}/ 2 f°°dp 



cos 



ir l / 2 El /2 

,3 

- - 2 7/3 p 







(F2) 



M") \ _ 1 (vqV?uT 

ISO I 64-27/12^ ^3/2^3/4 

" 32^ ( u_\ V2 



3%/3 \u h J 



(F3) 



h(x,y) 



in which the random magnetic field enters through u>h- 

3/2 

The argument of sine contains the term oc ui h , which 
can be presented as sho/h, where the constant s is equal 
to 16(2w/3wo) 3 ^ 2 - The factor in front of sine contains 
oc ft 3 / 2 . Then the gaussian averaging over h can be 
carried out analytically using the fact that s 1. This 
yields 



h 3/2 sh 



~h~ + 4 



hTs 1 ' 2 . f^' 2 s 2 ^ 



h(x,y) 



V6 



x exp 



^(^3-) 



3 S 2 / 3 

' 2 5/3 



(F4) 



Combining Eqs. (|F4|) and (|F3[) . one reproduces Eq. (|8.9|) 
of the main text. 
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